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ABSTRACT 


Understanding the movement of sediment in the nearshore region due to wave motion 
and longshore currents is important in beach erosion studies, and has tactical significance 
in beach front mine warfare. In the surf zone, air bubbles and sediment are both suspended 
within the water column. At the Naval Postgraduate School in Monterey, California, a 
sediment flux probe has been developed to study small scale processes. Using ultrasonic 
acoustic backscatter, the Coherent Acoustic Sediment Flux Probe (CASP) is capable of 
tracking the movement of scatterers within the surf zone. As it 1s important that the CASP 
system 1s capable of distinguishing between sediment and entrained air bubbles, laboratory 
experiments were run to determine the ultrasonic acoustic backscatter characteristics of surf 
zone bubbles. Bulk void fraction and optical sizing methods were explored to develop a 
means of measuring bubble populations produced in the laboratory for calibration of the 


backscattered energy received by the CASP system in the presence of bubbles. 
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I. BACKGROUND 


The natural movement of coastal sand has long been documented and studied by 
coastal engineers, oceanographers and developers. The ability to predict and forecast the 
movement of sand and other bottom sediments within the surf zone can be environmen- 
tally, economically and tactically significant. In the development of coastal communi- 
ties, man's efforts to protect and dredge channels and harbors has often adversely affected 
the longshore transportation of coastal sand and the natural evolution of barrier islands. 
Quite often, the efforts to protect man's investment in real estate and coastal structures 
has an adverse effect on the surrounding environment. During military operations, offen- 
sive and defensive mine warfare efforts can be seriously affected by the movement of 
sand. 

The sand budget in a given coastal region is governed by a combination of sinks 
and sources of sand. Sinks can be either natural (submarine canyons, dunes, and aeolian 
transport), or man made (sand mining and harbor improvements such as jetties). Sources 
are primarily natural, in the form of cliff erosion and river runoff. As previously stated, 
natural sources and sinks balance each other reasonably well until man's influence dis- 
torts the balance. 

Sand is transported throughout the surf zone by a variety of forcing mechanisms. 
Most dominant among these is movement due to wave forcing and longshore currents 
induced by deep water wave fronts breaking along the coastline obliquely. This forces a 
longshore component of momentum flux which, in turn, results in a longshore mass 
transport of sediment. 

Locally, sand is known to move on and off shore, in a seasonal cycle, on some 
coastal beaches. For example, Carmel Beach, in central California, reaches maximum 
width in September and October. However, by April, it 1s quite narrow, almost non- 
existent. The sand which makes up the scenic and famous expanse of the beach in the 
late summer and early fall is transported off shore and deposited in a sand bar. The sand 
bar seems to be relatively constant and predictable in its location and depth. It is 
reasonable to assume that this movement of the beach sand is dependent on the 


frequency and amplitude of the waves, which are significantly different in the winter and 
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summer. Tactically, this offshore movement and formation of bars can be very 
significant during military operations. In the area of mine warfare, this movement can 
negate offensive mining efforts by scouring the concealing sand from around buried 
mines. This can, of course, lead to much simpler hunting and sweeping efforts. The 
effort of planting mines on a beach to contain an enemy force could be wasted if the 
beach is likely to vanish completely in a matter of weeks. Defensively, it is quite 
obviously easier to hunt and sweep exposed mines than buried mines. On the other hand, 
proud mines could end up being buried by the deposit of ephemeral beach sand, requiring 
that the time of planting these mines be chosen appropriately. 

The underlying physical processes responsible for the movement of sand and 
sediment within the surf zone have been widely studied. However, the physical study of 
the movement of sediment within the surf zone is quite difficult. The surf zone itself is a 
highly turbulent and physically complicated region. A primary hurdle to scientific 
understanding of this region as it relates to sediment transport is the opacity, both visual 
and acoustic, of the water. The presence of bubbles entrained by breaking waves makes 
it difficult to track sediment movement. Without the presence of bubbles, it would be a 
relatively simple matter to study sand movement visually from above or, or acoustically 
in situ. in reality, however, the presence of bubbles makes visual tracking of sand nearly 
impossible and acoustic tracking quite difficult. Tactically, the opacity of the region 
makes it very difficult to detect and sweep mines within the surf zone. The presence of 
Strong scattering conditions due to high concentrations of both sediment and bubbles 
severely reduces the utility of acoustic mine countermeasures systems. 

There is a need to measure the physical processes contnbuting to surf zone 
sediment movement, which can resolve the small scale turbulent processes responsible 
for sediment suspension and transport. These studies would provide data concerning the 
nature of sand transport, and could eventually lead to tactical improvements in the areas 
of surf zone mine warfare and mine countermeasures. In an effort to develop a method 
to track and predict the movement of coastal sand, a series of experiments were 
conducted on the beaches of Duck, North Carolina, a few miles north of Kitty Hawk. 
The experiments of DUCK 94 utilized the Coherent Acoustic Sediment Probe (CASP) 
developed by Professor Timothy P. Stanton at the Naval Postgraduate School in 
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Monterey, California to conduct acoustic observations of the movement of entrained 
sediment in the surf zone, (Stanton [1993]). The CASP system uses ultrasonic acoustic 
backscatter information to track the movement of scatterers (sand particles, bubbles, and 
biotics in the water column). The capabilities of the CASP system will be discussed in 
detail. A difficulty with using ultrasonic sound in the surf zone is that both air bubbles 
and sand have the ability to scatter the sound energy, although with different 
characteristics. The objective of this thesis is to develop a discriminator to distinguish 
acoustic energy backscattered by entrained air bubbles from scattering by entrained 


sediments to measure sand and sediment fluxes at turbulent scales. 





Il, EXPERIMENTAL PLAN AND DEVELOPMENT 


In order to distinguish between acoustic backscatter from air bubbles and 
sediment, a method to measure the effects of bubbles on the Coherent Acoustic Sediment 
Flux Probe (CASP) backscatter level measurements was developed in a series of 
laboratory experiments. This required a means to produce a distribution of bubbles in a 
test volume, and a method to size and count the bubbles produced. The primary 
objective was to develop a method to discriminate bubbles from sediment, using dual 
frequency backscatter level and attenuation characteristics, and Optical Backscatterance 


Sensors (OBS) measurements. 


A. EXPERIMENTAL OUTLINE 

1. Development of Bubble Production System 

A bubble production system was required which was both simple to implement 
and capable of producing bubble clouds with size spectra similar to those in the ocean 
and of controllable densities. The production system had to produce a distnbution of 
bubbles in the test vessel whose total density could be easily controlled while having a 
size spectrum representative of ocean / surf zone conditions. Implementation was 
constrained by the physical size of the test vessel, a plexiglass tube 19 centimeters in 
diameter. A stirmng mechanism was used to create a homogeneous population of 
bubbles, and a conductivity cell directly measured void fraction within the test volume. 


The specific range of sizes required 1s discussed in a later section. 


2. Selection of a Bubble Measurement System 

As with the production system, a method of estimating the bubble size spectrum 
was required to be simple to implement, and to have a wide size measurement range and 
not perturb the bubble field. To allow robust size spectra to be estimated, an automated 
method of sizing was desired. Manual and visual methods of sizing and counting 
bubbles within the population were deemed to be too cumbersome for easy 


implementation. 


3. System Testing in Fresh Water and Salt Water Filled Tanks 
Controlled laboratory experiments were conducted in both fresh water and salt 
water. Initial tests were conducted in fresh, deionized water in an effort to reduce the 
influence of biologic impurities. This reduced the effectiveness of the conductivity cell, 
but provided a good baseline comparison with later tests using salt water. Sea water was 
not used, as control of biological contaminants was again a major concern. These tests 
were broken down into four major subcategories: 
a. Clear, bubble-free tank. 
b. Bubbles only. 
c. Sediment only. 


d. Bubble and sediment mix. 


B. EXPERIMENTAL SET UP AND DEVELOPMENT 


A literature search revealed many methods to both produce and measure air 
bubbles in a laboratory test tank. It was hoped that a proven method for both production 
and measurement could be utilized, to reduce the amount of time involved with the 
development of a laboratory experiment and apparatus, and mimimize the need for 
validation of either the production method or the chosen means of determining bubble 


population characteristics. 


1. Ocean Bubble Spectra Measurements 

The primary cnitena for selecting a bubble production system was that it simulate 
bubbles produced in the ocean during the spilling and breaking of waves, either in open 
ocean or in the surf zone. Several studies of the phenomenon of bubble production 
during the breaking of ocean waves have been performed. Published reports of field 
measurements of bubble populations in seawater have found a fairly consistent form for 
the entrained bubble spectra. Wu [1981] published a comprehensive summary of the 
results of several field measurements, most notably those of Medwin [1970], Kolovayev 
[1976] and Johnson and Cooke [1979]. Table (1) from Wu, shows a summary of field 


experimental conditions. 
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Table (1) Summary of Filed Experrmental Conditions (after Wu [1981 ]) 


Radius, um Radius, um 


Wu compared not only of the results obtained by these three investigators, but also of the 
limitations of each of the selected measurement systems. In order to obtain a meaningful 
comparison, Wu limited his study to size ranges common to all three data sets, 17 - 176 
um. The lower limit of this range was set by the results of Medwin and Johnson and 
Cooke, while the upper limit of this range, 176 tm, was imposed by Medwin's 
measurement techniques. Wu normalized the spectra measured by the three investigators 
to obtain a probability density function (pdf) of bubble radii. Figure (1a) shows the pdf 
of Kolovayev's results with a distinct peak near 80 um. Figure (1b) shows Johnson and 
Cooke's results with a peak in the spectrum near 50 um. Figure (1c) shows Medwin's 


results with no distinct peak between approximately 15 - 200 ym. The various plots 
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Figures (1a), (1b) and 1(c). Distribution of Measured Ocean Spectra (from Wu [1981}) 
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shown depict spectra measured at different sample depths. These figures clearly show 
the depth dependence of the number of bubbles present. In open ocean conditions, the 
entrainment level is shallow, although highly dependent upon wind speed and turbulent 
mixing. Deeper sampling levels result in fewer bubbles sampled, as shown in Figures 
(la) and (ib). Kolovayev measured bubble spectra at 1.5, 4 and 8 meters depth. Johnson 
and Cooke measured spectra at 0.7, 1.8 and 4 meters. Medwin conducted measurements 
at 1.5, 4.6, 7.6, 10.7, and 13.7 meters. Direct comparison of the pdf of Kolovayev's 
results with those of Johnson and Cooke, at 4.0 meters, shows that the peak of the 
distribution of the bubble radi for both investigators is near 80 um. Figures (2a) and 
(2b) show plots of bubble spectra at various depths and under different wind conditions 
as obtained by Kolovayev and Johnson and Cooke respectively. Figure (2c) shows Wu's 


plot of combined results without dependence on depth or wind velocity. 
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Figure (2). Ocean Bubble Spectra at Various Depths (from Wu [1981 ]) 


Closer examination of Figure (Ic) shows that Medwin's data does, in fact, have a plateau 
near 80 tm radius. While this does not represent a peak in his distribution, it does 


correlate with peaks shown in the data shown in Figure (2). The difference in magnitude 
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of the three plots can be attributed to errors associated with the sampling volume of the 
methods used to measure the bubbles and the smallest bubble resolved size. Medwin's 
method, using resonant acoustic scattering, has a much larger potential sample volume 
than the other methods. 

Su, et al, [1988] measured microbubbles in the North Sea optically, using dark 
field specular reflection, in a manner similar to the Light Scattering Bubble Counter tech- 
nique, Su et al [1994]. At a depth of 2 meters, they found microbubbles to fall within a 
range of 20-400 um. Figure (3) shows a comparison of results reported by Su et al 
[1988] with results reported by Johnson and Cooke [1979]. Although the measured 
concentrations do not compare too well below 50 pm radius, the spectrum shape and 
concentration levels compare well at large bubble sizes. 

Using an Acoustic Resonator Array, Medwin and Breitz [1989] and Breitz and 
Medwin [1989] reported results which were consistent with those of Su, et al [1988]. 
Figure (4) illustrates the close comparison between the results of Medwin and Breitz and 
Johnson and Cooke [1979]. Medwin and Breitz [1989] found bubble populations under 
spilling breakers at sea to have diameter of 25 - 270 um, but very significantly, no 
spectral roll off at small sizes. 

These field measurements all show distributions within the range of 
approximately 20 - 400 um despite the use of a variety of different measurement 
techniques, but with small radius cutoffs often influenced by the measurement technique. 
Details of the operation of the photographic and acoustic techniques, as well as others, 
will be discussed in detail later. 

Prior results of bubble generation in freshwater and seawater filled test tanks 
were researched to determine the diameter range which was attainable in the laboratory. 
Additionally, this would provide some insight into the effects of freshwater and seawater 
on bubble production, size and residence time. Glotov, et al [1962] determined a lower 
bound of bubble diameters to be 50 pm in a freshwater tank with mechanically generated 


wind waves. Bubble sizes were determined using a bubble trap, similar to that used by 
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Figure (3). Bubble Concentration Comparison (from Su, et al [1988]) 


10 


n(a) 
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Figure(4). Comparison of Bubble Spectra Under Breaking Waves 
(from Medwin and Bnietz [1989}). 


Kolovayev [1976]. In laboratory observations of spilling breakers in a controlled test 
tank, Medwin and Daniel [1990] measured bubbles between 50 um and 7400 pm 
acoustically. Cartmill and Su [1993], also using an acoustic resonator measured bubble 
size distributions in a test tank, ranging in size from 34 - 1200 um, under both freshwater 
and saltwater waves. 

Cartmill and Su [1993] found that bubble plumes generated in freshwater and 
saltwater were qualitatively different. Figure (5) shows that although the shape of the 
distribution curve 1s similar, bubble production 1n saltwater yields an order of magnitude 
more bubbles at a given radius. Saltwater bubbles were far more numerous at a given 
radius than those measured in freshwater. They found that this was primarily due to the 
easy coalescence of freshwater bubbles “whereas saltwater bubbles repel each other as a 
result of their surface physical-chemical properties". Scott [1975] found that bubbles in 
saltwater have a longer residence time due, in part, to "the mild surface activity of the 
salt itself. The result of the surface activity is to retard the coalescence of very small 
bubble produced by breaking waves, decreasing the velocity of mse of the entrapped air 
and increasing the time during which bubbles are observed to reach the surface and 
break." Haines and Johnson [1995] found that bubbles produced in saltwater were more 
numerous, smaller, and resident longer than in freshwater. It was therefore expected 


that, for a given bubble production rate, the bubble population in a saltwater filled tank 
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would have a higher concentration at small diameters and a lower concentration at larger 


diameters. 
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Figure (5). Bubble Size Distribution Comparison, Fresh and Saltwater 
(from Cartmill and Su [1993]) 


Based on these previously observed bubble size spectra, a size range of 20 - 400 
11m was determined to be a representative and desirable range for the bubble production 


system to be constructed. This size range would closely match the range of diameters 


found in field measurements at sea as well as those measured in laboratory experiments 


in wind wave tanks. 


2. Ocean Spectra Theory 

Bubbles are formed in the ocean primarily due to downward forces generated by 
breaking wind waves, or by shoaling and breaking waves in the nearshore. Ocean waves 
break for two reasons: kinematic and dynamic instability. Kinematic instability results 
when the velocity of individual particles of water within the wave becomes larger than 
the phase velocity, Cph, of the wave itself. Dynamic instability results when the 
acceleration of individual particles becomes greater than the downward force of gravity. 
In each case the particles separate from the leading edge of the wave and fall. The force 
of the particles falling as a group onto the leading slope of a wave, as shown in Figures 
(6a) and (6b), can entrain air bubbles below the surface of the ocean. These entrained 
bubble clouds remain beneath the surface after the wave travels along and are left to 
evolve. Koga [1982] stated that "the main mechanism of bubble formation in [wind] 
waves is intermittent bubble entrainment by an ordered convergent flow on the leading 
slope near the crest. Koga described this entrainment region as the area in which 
particles from the breaking portion of the wave impinge upon the leading slope of the 
wave. Provided that the downward force of these particles 1s sufficient to overcome the 
surface tension of the leading edge, bubbles will be entrained under the under the wave. 
With this conceptual model, the higher the wave, the greater the force, and the deeper the 
entrainment. If the area behind the wave entrainment zone 1s relatively quiescent, as in 
open ocean conditions, the bubble cloud will disperse relatively quickly as large bubble 
rise to the surface. If the area behind the breaking wave is highly turbulent, such as in 
the surf zone, the entrained cloud may be mixed through the water column several times 
before rising to the surface. This increased residence time can be expected to increase 


the size spectral levels of small bubbles. 
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Figure (6a). Breaking Wind Wave Bubble Entrainment Pattern 
(from Koga [1982]) 
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Figure (6b). Mean Distribution of Bubbles along a Representative Wave 
(from Koga [1982]) 
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3. Laboratory Bubble Production Methods 

Many methods of bubble production in a test volume were researched in the 
literature. Previous experimenters have produced streams using hypodermic needles fed 
by nitrogen tanks and electrolysis. Bowyer [1992] used measured volumes of water 
dropped from a height into the test volume. Bubble entrainment has also been 


accomplished through the use of surface and underwater jets. 


a. Electrolysis and Hypodermic Needles 

Bubble production via electrolysis typically uses a bubble - producing 
wire wrapped around a rod which sweeps through an arc at the bottom of the tank. This 
provides a vertically homogenous bubble plume, but bubbles produced at the ‘turning 
points’ in the arc of the arm are larger than those formed during the arm's sweep. This 
method was rejected as not producing a bubble population with significant size range, so 
it was likely that the bubbles created would not represent a statistically significant 
simulation of the spectrum estimated to be entrained in the surf zone. Additionally, 
since electrolysis produces hydrogen (H2) and oxygen (Oz) the bubbles will be highly 
soluble. The gas solubility effects the residence time of individual bubbles and can 
therefore effect the size distribution of the produced bubble population. Su, et al [1994] 
used a hypodermic system constructed of several needles mounted on a tube sweeping 
through an arc at the bottom of the test tank, very similar to the mechanism involved 
with the electrolysis method. Hypodermic needles do not produce the highly soluble 
bubbles, but tends to produce a narrow range of source bubble sizes, and was rejected for 


this reason. 


b. Underwater Jet Air Entrainment 

This method utilizes an under water nozzle acting as an eductor which 
sucks air into a water stream due to the velocity of the water within the nozzle. The 
entrained air is fed either from a compressed air supply or drawn from a gas reservoir. 
The size of produced bubbles varies from tens of microns to several millimeters. 
Authors found that this gave a more continuous bubble density over the range of 30 - 


1200 xm. This range of radii was attractive as it nearly met the criteria for matching in 
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situ measurements of bubble spectra under breaking waves [Breitz and Medwin, 1989]. 
Hypodermic, Electrolytic and Underwater Jet methods were briefly described by Su, et al 


(1994) and were used in their comparison of various bubble measurement systems. 


c. Above Water Jet Air Entrainment 

Open ocean tests have found fairly consistent results showing bubble 
spectra between 20 - 400 tum [Breitz and Medwin, 1989], as shown in Figure (4). A 
production method was required which would produce a continuous density, broad 
spectrum of bubbles which could be mixed to create an homogeneous bubble population. 
In order to do this, a bubble production method developed by Koga [1982] was chosen. 
This method involved an above water jet impinging upon the surface of the test tank. 
The jet entrains air bubbles due to the downward force acting upon the surface of the 
fluid in the test volume. Bubble entrainment is dependent upon overcoming the surface 
tension of the fluid. Koga {1982] found that this method, for his experimental 
parameters, created air bubbles ranging in size from 25 - 4500 um in diameter. This 
method is very similar to the under water jet entrainment method noted previously, but 
does not require an eduction nozzle. It also reduced the complexity and amount of 
tubing required in our experimental apparatus. Given the nature of set up and the lack of 
physical space within the test vessel to be used in the laboratory, this was an :mportant 
consideration. However, production of a bubble cloud is highly dependent upon the 
angle of impingement of the jet upon the surface, as shown in Figure (7). The figure 
clearly shows that at shallow incidence angles, as in Figure (7a), the entrained bubble 
cloud contains many relatively large bubbles, and the entire cloud remains near the 
surface of the tank. At steep angles of incidence, such as in Figure (7d), very few 
bubbles are created. Koga found that at moderate angles, the entrained bubble cloud had 
a broad spectrum of bubble radii, and the produced cloud was allowed to evolve farther 
below the surface than with shallow angles, which improved mixing within the stirred 
test vessel. The critical impingement angle was a very real concern, however. Due to 


the small size of the test vessel, flexibility in selecting an impingement angle would be 
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constrained. This was overcome through the use of a fairly short ngid nozzle fed by 
flexible Tygon tubing. This provided the necessary flexibility to select an optimum 
impingement angle without restricting the use of other apparatus at the surface of the 


tank due to the size of the nozzle apparatus. 
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Figure (7). Bubble Entrainment for Various Nozzle Inclinations (from Koga [1982]) 
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4. Bubble Production System Development 


Initial development and design of an above water jet for the production of air 
bubbles in the test volume closely followed the design used by Koga [1982]. A self 
priming water pump was initially used to circulate water from the test vessel, through a 
nozzle and back into the test vessel. Imitial tests were conducted by visual inspection of 
the produced bubble cloud. 

Visual testing of the water jet showed that this particular method was both easy to 
set up and easy to control. Varying the size of the nozzle appeared to alter the spectrum 
of the bubbles produced by varying the volume of the flow, while adjusting the flow 
velocity of the jet for a given nozzle diameter with a needle valve, altered the total 
density of the produced cloud without noticeable affect on the produced spectrum. 
Initial tests of this system, without direct measurement of the resulting bubble 


population, used three nozzle diameters as shown in Table (2). 


Nozzle Diameter 
.046 in (.117 cm) 
.082 in (.208 cm) 
.183 in (.465 cm) 


Table (2). Entrainment Jet Nozzle Sizes for Initial Tests 


During tests it was found that the produced spectrum was both dependent on 
nozzle diameter, and therefore flow velocity, as well as the angle of impingement. These 
results are consistent with those reported by Koga [1981]. He found that for an 
impingement angle, 0, of 45 degrees, the distribution of the produced spectrum was 
broader than at 0 = 30 degrees and 0 = 50 degrees. Both of these angles produced a 
narrower distribution of sizes. In addition to a dependence on 9, bubble production was 


also highly dependent on the height of the jet above the surface of the water. As reported 
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by Koga [1981], the greater the distance above the water, the more unstable the flow 
from the jet became and the more confused and unstable the produced bubble population 
became. It was found that to produce a stable, easily controlled bubble population, the 
nozzle could be positioned approximately 1 centimeter above the surface of the water. 
Above this height, the produced cloud contained a great many very small bubbles, but not 
a large range of bubble sizes. If positioned too near the surface, the nozzle had a 
tendency to become submerged due to surface agitation brought on by the laboratory 
stirring apparatus used to produce a homogeneous bubble cloud within the test vessel. 
The pump used in the initial tests was a March Manufacturing model MDX 115 volt, 
1/55 horsepower pump. 

Due to the low flow rate of the pump, it was necessary to reduce the interior 
diameter of the tubing to produce a flow velocity capable of overcoming the surface 
tension of the water in the test volume. The exit diameter of the pump was 0.5 inches 
(1.27 centimeters). This was decreased at the nozzle to the diameters shown in Table 
(2). It was found that the 0.183 inch diameter nozzle produced a low velocity jet which 
was only capable of producing a narrow spectrum of fairly small bubbles. The 0.082 
inch diameter nozzle had a broader spectrum, but produced a very low density of 
bubbles. Since these initial tests were conducted with the velocity controlling needle 
valve fully open, it was essential to produce a high density, broad spectrum population in 
order to have a workable density population at lower flow rates. The final nozzle, 0.046 
inches produced a broad spectrum with a high density population. Qualitative 
observations of the produced bubble clouds revealed that this method produced clouds 
which had controllable density with a relatively constant size spectrum. The critical 
requirement was to produce broad spectrum clouds which had controllable density but 
which maintained a constant spectral range, unaffected by flow velocity. A diagram of 


the production system and associated apparatus is shown in Figure (8). 
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Figure (8). Above Water Jet Entrainment Apparatus 


5. Bubble Size Spectra Measurement 


Several criteria needed to be met in selecting a method of bubble size spectra 
measurement. The measurement method had to be capable of measuring bubbles 
without disturbing the bubble cloud. This would restrict the use of any method which 
physically trapped a single bubble, or disrupted the natural evolution and rise of any 
individual bubble within the population. The sizes of the measured bubbles should not 
be restricted by the physical limitations of the selected equipment. Perhaps most 
important among the considerations were those of cost and ease of implementation. The 
system had to be relatively inexpensive and simple to operate and install. Several 
methods of bubble measurement were investigated in a literature search, including: 
Light Scattering Bubble Counters (LSBC), Photographic Bubble Imaging System (PBIS), 
Acoustic Resonator Arrays (ARA), and Laser Holography and Coulter Counters. LSBC, 
PBIS and ARA systems were compared for use In sizing bubbles by Su, et al (1994). The 
goal was to select a measurement method which was robust, easily implemented, had a 
large sampling volume, had the ability to measure bubbles of 20 - 500 um radius 
accurately, quick response time and could sample the bubble cloud accurately using 


automatic data analysis. Additionally, the method would preferably have the ability to be 
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deployed in the field to obtain in situ measurements. A review of methods described in 


the literature follows. 


a. Light Scattering Bubble Counter 

The Light Scattering Bubble Counter system uses dark field specular 
reflection for detection of microbubbles. White light is reflected by the air - gas 
interface of the bubble with the surrounding water. Typically, bubbles are drawn into a 
sampling volume by a pump, and reflected light intensity 1s detected by a photomultiplier 
at a predetermined scattering angle with respect to the incident light. Figure (9) shows a 
diagram of the apparatus used by Su, et al [1988]. With this particular set-up, Su, et al 
[1988] chose a scattering angle of 125 degrees in order to detect a strong forward 


scattering portion of the light energy scattered by particles within the sample volume. 
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Figure (9). Light Scattering Bubble Counter (from Su, et al [1988]) 
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The intensity of the reflected light, as detected by the photomultiplier, is proportional to 
the square of the bubble's radius. Su, et al [1994] calibrated his system by determining 
the size of the individual test bubbles from their terminal rise velocity. This system, as 
reviewed by Su, et al [1994], was limited by it's inability to sample bubbles of greater 
than 200 pm radius due to the physical size of the intake of the sampling volume. This 
limitation could easily be alleviated by changing the size of the sample volume entry. 
This system was rejected for a vanety of reasons. The limited size range of the bubbles 
measured, 10 - 150 um, and the small sampling volume were problematic, but its 
inability to sample more than a single bubble simultaneously was seen as severely 
limiting in its application to the high bubble densities to be generated in the laboratory 
tank measurements. There was also a concer that the method of pump sampling 
bubbles into the sampling volume might distort not only individual bubbles, but the 
overall spectrum by breaking up large bubbles in order to receive them in the sampling 
volume. Ling and Pao [1988] used a similar apparatus to measure micro bubbles in the 
North Sea. With a selected scattering angle of 125 degrees, they detected a pattern of 
specularly reflected light on the surface of air bubbles which they called “hot spots”. 
They determined that the size of this reflected light spot was directly dependent upon the 
size of the bubble and was related to the selected scattering angle. Figure (10) shows the 
geometric optical scattering pattern as seen by the photomultipliers mounted on Ling and 
Pao's equipment. 
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Figure (10). Specular Reflection off the Surface of an Air Bubble in Water 
(from Ling and Pao [1988]) 
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b. Photographic Bubble Imaging System 

Su, et al [1994] used a Photographic Imaging System in a companison of 
several measurement techniques. This system uses photographic images of bubbles 
illuminated by a combination of three strobes, separated by 120 degrees. This provides 
three bright points around the circumference of each individual bubble. Figures (11) and 
(12) show the orientation of the camera and strobes and the pattern of strobe reflections 
upon the bubble's surface respectively. In each image of a bubble, the three dots occur 
due to specular reflection of the photographic strobes. The sample volume was 


determined by the width and length of the image, and the depth of field of the lens. 





Strobe 
Figure (11). Photographic Bubble Imaging System (from Su, et al [1994}]) 
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Figure (12). Strobe Pattern on a Bubble's Surface (from Su, et al [1994]) 
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Manual counting and measuring of bubbles is performed with a microscope. Although 
this system is capable of sizing bubbles of up to 500 um, a lower bound of 50 pm is 
imposed by the finite resolution of the microscope, the sensitivity of the photographic 
system and the graininess of the film used in the photographic process. The upper bound 
is due to the low statistical probability of encountering a bubble larger than 500 um 
within a sampling volume, and the limited number of samples due to the laborious visual 
analysis of the photographs. This spectral range was limited, as expected, as much 
smaller radi (~10 wm) bubbles would be encountered both during laboratory 
experiments and during in situ measurements. However, the lower bound could have 
been improved by using a more powerful microscope. Results obtained by Su, et al 
[1994], using this system in the laboratory, correlate with the in situ findings of Walsh 
and Mulhearn (1987). This system does not seem to disturb the flow of the bubbles, but 
Su, et al found it only practical for calibrating other systems in the lab due to the 
extremely labor intensive manual microscopic measurements of the photographic 


negatives. 


c. Acoustic Resonator Array 

The Acoustic Resonator Array provides a continuous distribution of 
bubble sizes without tedious manual counting and microscopic measurements. The 
system developed by Breitz and Medwin [1989] consisted of a transducer and reflector 
plate separated by 20 cm. When a broadband acoustic source excites the volume 
between the plates, standing waves arise due to a series of harmonic resonances governed 


by the relationship: 


f=enc/2L 
Equation (1) 


where c = speed of sound, n = number of the acoustic harmonic and L = distance 


between the plates. This results in a cavity fundamental frequency of approximately 3 


IS, 


kHz, dependent upon the speed of sound in the water, with multiple harmonics over the 
bandwidth of the white noise transmitted into the cavity. This yields 1 fundamental and 
32 harmonics over a 100 kHz band. Resonant bubble frequency is given by: 


fo=(1/a)¥(3yPo/eo) 


Equation (2) 


where a = bubble radius, y = ratio of heat capacities, P, = hydrostatic pressure and po = 
the density of the water. Over the 100 kHz white noise range, this relationship covers 
bubble diameters from 1200 4m down to 34 um, corresponding to resonant frequencies 
between 3 and 99 kHz. Higher frequencies could be used, but the extinction cross 
section, 6, 1s no longer due simply to the resonance of the individual bubbles but to the 
increasing viscous effects in smaller bubbles. The separation between the plates could 
also be altered to expand the range, but then the resonance quality factor, q, 1s lowered. 


The bubble size distribution is given by the formula 


n.(a) = (2 x Af) / (co 10° 8a) 
Equation (3) 


where n, = bubble density, c = speed of sound, and Af = the change in resonance width 
calculated from the sound frequency spectrum comparison with bubbles present and 
absent. A similar system, as described by Su, et al [1994], is capable of sampling bubble 
radii of 34 - 1200 ym. It also samples the bubble cloud with a volume of 1250 cc, a 
much larger sampling volume than any of the other systems investigated. Both of these 
factors, coupled with the ease of use both in the field and the laboratory made this highly 
attractive. Unfortunately, the equipment necessary to employ this method was 
unavailable for the laboratory portion of the experiments, and the small test vessel 


precluded using borrowed existing field systems. 
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d. Laser Holography 

O'Hern, et al [1988], used a pulsed ruby laser directed through a sampling 
volume to obtain holographic images of scatterers within the volume. A portion of the 
beam is diffracted by particles, while the remainder passes through the volume. High 
resolution film records the interference pattern of the diffracted and undiffracted laser 
light, and the image is subsequently analyzed to size particles illuminated within the 
sampling volume. The developed hologram is illuminated with a He-Ne laser beam. 
This image provides a three dimensional representation of the illuminated bubble cloud. 
Particles are discriminated from bubbles by visual observation. Bubbles are brighter, 
more rounded and have a central dark spot due to internal focusing of the laser light. 
This central dark spot is due to the fact that O'Hern viewed the illuminated bubbles along 
the axis of the beam. Due to the design of our laboratory tank, and the many 
complexities of laser holography, this method proved to be impractical. O'Hern's laser 
holography system has extremely good object resolution, approximately 5 um. Again, 
the analysis of produced images would be very tedious as objects within the sample 
volume would be sized by visual measurement only, rather than an electronic, or 
automated means. During ocean tests of this system, O'Hern, et al, found holographic 
techniques to be a reliable method of bubble measurement, but surprisingly no bubbles 


larger than 50 microns were detected. 


e. Coulter Counters 

Coulter Counter systems detect changes in conductivity as particles, or 
bubbles, in a fluid pass through a detection volume. The amplitude of the signal is nearly 
proportional to the volume of the particle, as long as the particle dimensions are less than 
40% of the orifice diameter. O'Hern, et al [1988] used a Coulter Counter apparatus 
with an orifice diameter of 140 um, limiting the maximum diameter of measured bubbles 
to approximately 50 pm. It is, of course possible to use a larger orifice diameter to 
enable larger bubbles to enter the counter. A Multi-Channel Analyzer sized incoming 


signals and sorted them into 128 separate channels. It was noted by O'Her that Coulter 
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Counters may suffer an order of magnitude error in estimation in particulate content. 
This is possibly due to the disintegration of the particles as fluid is drawn into the 
measurement volume, therefore, the measured spectra of bubbles would not necessarily 
be accurately represented. Large bubbles could easily be broken up as they passed 


through the orifice to be counted, making it an unattractive option. 


6. Measurement Apparatus 

Following this literature review and careful consideration, a bubble measurement 
system was decided upon which consisted of a hybrid of several of the previously noted 
systems. The measurement system consisted of a high intensity laser beam which was to 
be used to illuminate individual bubbles within the bubble cloud, a microscope and 
video camera assembly, and framegrabbing software. This system utilized the bubble 
illumination theories of the Light Scattering Bubble Counter (LSBC) system. By 
capturing a video image, Photographic Bubble Imaging System (PBIS) fundamentals are 
exploited, but analysis was enhanced through the use of automated image processing 
software. The narrow beam, high intensity of a laser was utilized without the 
cumbersome nature of holographic imagery and the complexity of holographic apparatus. 
The bubble population would not be disturbed and would evolve without physical 
constraints as noted in research of bubble traps, Light Scattering Bubble Counters and 
Coulter Counters. Lasers have been used in many instances to measure bubbles and 
other suspended objects in laboratory and field experiments as previously noted. In most 
cases, however, lasers have been used to measure single objects rather than large groups 
or clouds of objects. 

Initially, a beam spreader was mounted on the head of the laser to create a fan of 
light spread over 30 degrees, with a depth of 0.3 mm. This fan of light was to be 
reflected by a mirror apparatus as shown in Figure (13). It was thought that this 
geometry would provide bubble images similar to those obtained using a Photographic 
Bubble Imaging System and as shown in Figure (12). Only bubbles within the fan of 


light would be illuminated, but none outside of it. This provided a well defined sample 
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light would be illuminated, but none outside of it. This provided a well defined sample 
volume and high illumination intensity for the video camera and microscope. The laser 
was aimed down through the surface of the test vessel with the camera positioned to view 
the illuminated cloud normal to the plane of illumination. However, it was found that 
the surface disturbance created by the stirring action within the test vessel coupled with 
the reduced light intensity caused by the beam spreader made this set up ineffective. 
Additionally, within the turbulent, high density bubble cloud within the test vessel, there 
were many instances where each bubble did not have three distinct specular "hot spots” 
in the captured image. This was primarily due to the interference of other bubbles 
blocking the light rays from stnking the surface of another bubble. Therefore, bubbles 
recorded by the video camera had between one and three "hot spots". The random nature 
of the representation of a bubble by "hot spots" resulted in this optical geometry being 
rejected. Instead, a single beam of light was projected into the side of the tank without 
the attachment of either a beam spreader or the mirrors. This prevented the surface 
interference, light intensity and random “hot spot" problems without eliminating the 


ability to attain a 90 degree scattering angle within the test vessel. 
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Figure (13). Laser and Mirror Assembly 
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7. Optics Theory 


A simple model of the reflection of an incident plane of light on an air bubble 
suspended in a volume of water on a ray of incident light was developed to understand 
the imaged reflection measured by the laser / microscope system described in section 6. 
When a light ray strikes the surface of a spherical, clear, or transparent object, some of 
the energy 1s reflected at an angle equivalent to the incident angle, t. In this case, t = 90 
denotes central incidence and t = 0 denotes grazing incidence. Figure (14), from Van de 


Hulst [1957], illustrates the effect of a single ray of light incident upon the surface of an 





Figure (14). Single Light Ray Trace for an Air Bubble in Water 
(from Van de Hulst [1957]) 


air bubble in water. On a sphere, where the slope of the surface becomes increasingly 
steep as One travels outward from the center, this means that parallel light rays striking 
the surface are scattered over a wide range, as shown in Figure (15). This figure shows a 
two dimensional slice of an air bubble in water. The incident light rays have been 
generated by a distant source, and all rays are assumed to be parallel. Most of the 
incident energy is reflected, as shown in Table (4). The remainder of the light energy 1s 
transmitted into the bubble and refracted prior to striking the interior surface. A portion 


of this energy is, in turn, internally reflected, and a portion is transmitted out of the 


i} 
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bubble, refracted again by the surrounding fluid. The energy departs the bubbles surface 


at an angle equal to the incident angle, t. Energy refracted in the interior of the 
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Figure (15). Multiple Light Ray Trace for an Air Bubble in Water 


bubble follows Snell's Law: 


cos t'/n,=cos t/n, 


Equation (4) 


where Tt’ is the internal refraction angle, n, and n, are the indices of refraction of water 


and air respectively. Van de Hulst gives this relationship as: 


cos t =(1/m) cos t 


Equation (5) 


where m is the ratio of the indices of refraction, m=n,/n, For an air bubble in water, 
Van de Hulst gives the index ratio, m, the value 0.75. Halliday, et al (1993), give the 


index of refraction for air and water as follows: 


Index of Refraction, n 
Air (STP) 1.00029 
Water (20 degrees C) 


Table (3). Indexes of Refraction for Light in Air and Fresh Water 
(after Halliday, et al [1993]) 










These values result in a value of m equal to 0.752. This physical relationship for internal 
refraction holds until the angle of incidence becomes equal to or greater than the critical 
angle for the given index of refraction. The critical angle, 0 crt, as defined by Halliday, et 
al [1993], is the angle of incidence above which total internal reflection occurs, i.e. for 


angles of incidence greater than 6 crit, rays will experience no refraction, only reflection. 


Critical angle is given by the relationship 


8 crit = arcsin( n, / n,) = arcsin(m) 


Equation (6) 


For a bubble of given radius, r, assumed to be a perfect sphere, the angle of incidence for 
a beam of parallel light rays is dependent on the slope of the exterior surface of the 
sphere. If one considers the case of a vertical light ray, normally incident on the bubble's 
surface at the center of the bubble, the critical angle occurs at 0.75 (r). Light rays 
emitted from a laser without divergent optics are all considered to be parallel. Given 
that the ratio of the indices of refraction, m, has a value of 0.75 [Van de Hulst, 1957], the 
critical angle is 48.59 degrees. A sphere of radius, r, cut directly through the center by a 


plane forms a circle of radius, r. For this circle, described by the equation 


xt+yap? 


Equation (7) 


OT 


yovir?23 


Equation (8) 


the slope of the surface at any point displaced from the center of the bubble by a 
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horizontal distance, x, given as a function of x as 


® . xe (Vfr?-x7)) 1 
dx 


Equation (9) 


For a vertical ray of light, the angle of incidence, t, at a distance 0.75 (r), measured from 
the center of the bubble, is exactly 48.59 degrees. Therefore, any light energy incident 
upon the surface of the bubble from 0.75 (r) to the edge will be completely reflected, as 
shown in Figure (15). Figure (15) shows a bubble viewed from the side, that is, viewed 
from an angle normal to the viewing angle of the microscope / video camera assembly. 
The parallel ray source producing the light rays 1s above the bubble, and all rays incident 
upon the surface of the bubble are vertical and parallel in this viewing plane. Since the 
camera will only see light scattered toward it, only the front half of the bubble 1s shown 
to be scattering light energy. In reality, the bubble would scatter energy from every point 
on its upper surface in addition to refracting energy through the entire lower half. This 
figure does not account for rays which are internally reflected one or more times and then 
exit the bubble toward the camera, since rays experiencing multiple interior reflections 


and refractions are strongly attenuated. 
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Table (4). Light Energy Distribution for an Air Bubble in Water 
(after Van de Hulst [1957]) 


Total reflection, no refraction 43.8% 


Table (4), above, clearly shows that most of the light energy incident upon the surface to 
a bubble is either specularly reflected or refracted twice, exiting the bubble as shown in 
Figure (15). These rays account for the distribution of 93.5% of the incident energy. For 
this reason, other rays, those experiencing interna! reflection, or more than two 
refractions, were neglected in preliminary optical ray tracing efforts. With the degree of 
magnification required, the amount of reflected and refracted energy becomes a concern 
due to the sensitivity of the camera. The addition of a video microscope also affects the 
amount of light required to be received by the camera. The shutter speeds required to 
freeze bubble movement in the image (between 1/1000 th and 1/10,000 th seconds), also 
require much brighter images than an unshuttered camera. As Figure (15) shows, light 1s 
scattered in many different directions as rays strike the surface of the bubble. Depending 
on the width of the viewing field of the camera, however, only a few of the rays are 
actually raged by the camera. Figure (16) shows the predicted pattern of total specular 
reflection from a sphere. It should be noted that this figure does not show any energy 
which is transmitted through the interior of the bubble, but only the specular reflection of 


light energy off the upper surface of an air bubble. 
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Figure (16). Specular Reflection of Light Energy Incident on a Bubble 
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Energy refracted twice through the intenor of the bubble would not be detected by the 
camera, as none of it 1s predicted to scatter in the direction of the camera. 

Given the theoretical scattering patterns of the light rays, it was predicted that the 
image produced of the illuminated bubble with a single illumination source would be a 
bright dot representing an area of specular reflection off the top surface of the bubble 
only. In order for the automatic sizing and counting software to detect bubble images, it 
would be necessary create an image of a bnght dot which was easily distinguished from 


the background. 


8. Geometrical Optical Approximations 

Air bubbles suspended in water have been measured by a variety of methods, as 
previously discussed. The majority of the methods rely on physical measurement of the 
actual bubble, either via void fraction volume, photographic image or mse velocity 
measurements, to name a few. Measurements using scattered light have often relied on 
illuminating a single bubble with strobe lights, such as in the Light Scattering Bubble 
Counter [Su et al, 1994]. All of these methods are only applicable for single bubbles, 
and would not be appropmnite in a high density cloud of bubbles typical of the surf zone, 
or which could be simulated under laboratory conditions. The use of light scattering 
method under these conditions requires a technique capable of sizing individual bubbles 
within a population. For this reason, a single direction illumination source was used, and 
the size of the specular reflection spot was determined from the captured video image. In 
order for the theory of geometrical optical approximation to apply, as described by Van 


de Hulst [1957] and Baldy and Bourguel [1985], the following restrictions are required 
d 
m(—)»1 
( 7 


Equation (10) 
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where d is the diameter of the air bubble and A = the wavelength of incident light, and, 


Aa210° 


Equation (11) 


where Aq is the solid angle subtended by the optic aperture of the viewing apparatus. In 
the size range prescribed by the first assumption, bubbles are characterized by a scattered 
light "hot spot" which is smaller than the bubble. The size of the "hot spot" is directly 
proportional to the size of the bubble (Baldy and Bourguel [1985] and Ling and Pao 
[1988]) 


Os CA bble 


Equation (12) 


where dhotspot 18 the diameter of the scattered light "hot spot", dbubbie is the diameter of the 
bubble and cs ts a constant of proportionality which 1s a function of the scattering angle, 
9, and the radius of the bubble. If the scattering angle is equal to the angle of incidence 
of the light energy, the "hot spot" will appear as a circle. The larger the scattering angle, 
the more elliptical the shape of the "hot spot" will become. Ling and Pao [1988] discuss 
the shape of the light source as seen by the detector. In their case, a 125 degree 
scattering angle was used between the axis of a white light source and a photomultiplier. 
They found that the diameter of the “hot spot" on the surface was approximately 0.10 
times the diameter of the bubble. Therefore, using previous notation, cs = 0.10. This 
value 1s somewhat larger than that found by Baldy and Bourguel [1985]. They utilized a 
laser light source with a 90 degree scattering angle and found cs = 0.06. Figure (10), 
from Ling and Pao [1988], shows the relationship between the size of the "hot spot” and 
the scattering angle. The upper portion of the figure shows a ray trace of the ray 


scattered at the critical angle and the ray scattered at the selected scattering angle. The 
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lower portion of the figure shows a representation of the image of the "hot spot" as seen 
by the observer. This portion of the figure has been rotated in order to simplify the 
picture. 

Figure (10) shows that if a shallower scattering angle, 90 degrees for example, 
was selected, the image of the "hot spot" would have a much smaller minor axis than the 
image shown. The more acute the angle between the critically scattered ray and the ray 
scattered at the selected scattering angle, the smaller the displayed “hot spot" image. 
This is the primary reason for the difference in the proportionality constant found by Ling 
and Pao [1988] and Baldy and Bourguel [1985]. Baldy and Bourguel [1985] state that 
the aperture of the viewing apparatus used was a = 14 degrees. Although Ling and Pao 
[1988] do not detail the aperture of their photomultiplier, it is safe to assume that it 1s 
sufficiently large for them to be able to use the geometrical optical approximation to 
size bubbles in their experiment. 

A detailed description of the methods of video capture and data analysis follows 
in Chapter V. 
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Il. COHERENT ACOUSTIC SEDIMENT FLUX PROBE (CASP) 


The Coherent Acoustic Sediment Probe (CASP) is a multi-frequency underwater 
sediment flux sensor capable of measuring both the turbulent scale velocity field and 
suspended sediment flux. It was initially designed and developed at the Naval Postgraduate 
School (Stanton, 1993) to conduct studies of turbulence in the oceanic boundary layers. The 
head of the sensor contains a single 1.4 MHz transducer surrounded by three 5.3 MHz 
transducers, as shown in Figure (17). The 5.2 MHz transducers are positioned radially 
around the center transducer, separated by 120 degrees. Two axis tilt sensors and three axis 
accelerometers monitor sensor orientation and motion. Electronic processing and control 
modules transmit a high speed data stream to a shore based data collection station. 

The acoustic package of the CASP is capable of operating in either a monostatic or 
bistatic mode. In the bistatic mode, a single 5.2 MHz transducer emits a series of 32 short 
pulses, while the other two transducers receive backscattered acoustic energy from the water 
column. Following sampling of all thirty two pulses by the other two transducers, the 
transmission of the next sequence of pulses is electronically switched to a different 
transducer head. This sequence is repeated with each of the 5.2 MHz transducers. In the 
monostatic mode, a series of three somewhat longer pulses, referred to as a triplet, is 
transmitted from each 5.2 MHz transducer. The transducers switch between transmit and 
receive mode sequentially to measure backscattered acoustic energy amplitudes in each of 
a series of 1.68 cm range bins out to a range of 1.2 meters. 

For calibration purposes, the CASP was positioned in a 600 gallon acrylic test tank 
in the laboratory. The test vessel in which bubbles or sediment were contained was 
positioned such that the acoustic beams transmitted by the CASP head would project 
through an acoustically transparent window in the side of the vessel. This allowed 


backscattered energy from within the test vessel to be measured by the transducer head 


4] 


without nearfield interference problems. Figure (18) shows a schematic of the CASP 


orientation within the test tank in relation to the test vessel. 
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Figure (17). Schematic of CASP System Acoustic Beam Geometry 
(from Stanton [1993]) 
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1. COHERENT ACOUSTIC SEDIMENT PROBE (CASP) 
2. ACOUSTICALLY TRANSPARENT MEMBRANE 

3. ACRYLIC NOZZLE AND DUCTED PROPELLER 

4. ACRYLIC TEST VESSEL 

5. VARIABLE SPEED MOTOR 

6. OPTICAL BACKSCATTER SENSOR 

7. ACRYLIC LABORATORY. ACOUSTIC TEST TANK 


8. ENTRAINED BUBBLE CLOUD 


Figure (18). Schematic of CASP System and Test Vessel (from Stanton pewanh) 





IV. ACOUSTIC EFFECTS OF AIR BUBBLES IN WATER 


A. BACKGROUND 

Gaseous bubbles occur naturally in the ocean as clouds of air bubbles entrained by breaking 
surf or falling rain, gas filled bladders in fish and as lung cavities in marine mammals. In general, 
naturally produced air bubbles occur near the surface when produced by falling rain or breaking 
waves. The larger of these bubbles rise quickly to the surface, while the smaller bubbles tend to stay 
resident in the water column for long periods of time, particularly when turbulent mixing is 
occurring. The presence of high density clouds of bubbles, and therefore a significant amount of 
air, within a volume can have a profound effect on the speed of sound, and the attenuation and 
scattering of sound energy due to scattering within that volume. This effect is dependent upon the 
size distribution and spectral density of the bubble cloud. These influences will be examined in later 


sections. 


B. SOUND VELOCITY 

Air which ts dissolved in a volume of water has very little effect upon the speed of sound 
within that volume, (Unck {1983]). However, air which is present as suspended bubbles has a 
significant effect upon the speed of sound. Urick states that when the bubbles are much smaller 
than resonant size, sound speed is dependent upon the density and compressibility of the air and 
water and the void fraction present within the volume. Void fraction is defined as the amount of 
air, by volume, present within a sampled volume of water. The simplified relationship for velocity 


of sound within the mixture is given by Urick as: 


ve) L- ple! od. 
px "Nl 7+2.4x104P 


Equation (13) 
where vw = sound speed in water, p = density of the air water mixture, kK = compressibility of the 
mixture and B = void fraction of air. Even at small void fractions (0.01 percent), the sound speed 


can be reduced by as much as 53%. At high frequencies, above the resonant frequency of the 
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smallest bubbles in the cloud, however, the effect of suspended air bubbles on the speed of sound 
is negligible. Unck gives the relationship between the speed of sound in the mixture, c, to the speed 


of sound in water alone, Co, as: 


Po -2 47 na 
Co =6 
(2nfy 


Equation (14) 
where n = the number of bubbles in a unit volume and a = the radius of the bubble. Bubbles can be 
considered to be acoustically "large" compared to the resonant size, when all bubbles present within 
the sampled volume are much larger than the resonant bubble size. This is the case for the surf zone, 
given previously published studies of bubbles formed due to breaking waves, and the transmitted 
acoustic frequencies used by the Coherent Acoustic Sediment Flux Probe (CASP). Figure (19) 


shows a plot of expected surf zone bubble sizes and their corresponding resonant frequencies. 


C. SOUND IN THE OCEAN 

Sound propagation in the ocean is effected by a vanety of phenomena. As the Coherent 
Acoustic Sediment Flux Probe (CASP) works over ranges of 1-2 meters at 1.3 and 5.2 MHz, 
temperature, salinity and pressure can be considered constant and their effects on the speed of sound 
can be considered negligible. Regardless of range, however, sound energy can be attenuated due 
to spreading, absorption and scattering within the sample volume. Spreading of sound energy within 
the water column is considered to be spherical unless boundaries are encountered. Spherical 


spreading is governed by the relationship 
TL=20logr 


Equation (15) 
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Figure (19). Surf Zone Bubble Resonant Frequencies 
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where TL is defined as the Transmission Loss experienced by the sound energy due to spreading 
alone as it travels through the water, and is measured in decibels referenced to 1 micropascal (db 
re 1 Pa) at a range of one meter from the source on the acoustic axis. When losses due to 


absorption are included in the spherical spreading law, the relationship becomes 


TL=20logr+ar 

Equation (16) 
where a = the absorption coefficient in decibels / meter. The absorption coefficient can be 
determined using Thorpe's equation 


O.1f? 40f’ 


+2.75x10 ~*f? +0.003)x10° 
1+f? 4100+? 





a. =( 
Equation (17) 


where the frequency, f, is in kHz. Sound wave energy attenuation due to the presence of bubbles 


can be calculated with the relationship 


I, 
a ,=10log— 
fl, 


Equation (18) 


where I, and I, are reference intensities measured on the acoustic beam axis one meter apart. If x 


resonant bubbles are present within a volume, each with an extinction cross section, o,, the 
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attenuation coefficient becomes (The Physics of Sound in the Sea, 1947) 


aw ,=4.34n0, 


Equation (19) 


where the coefficient has the units of decibels per meter. However, since it is assumed that the 
CASP system will not excite resonance in any of the bubbles found to be resident within the surf 
zone, this entire term can be disregarded. Given that the CASP is essentially an active sonar emitter, 


the transmission loss due to range and weak attenuation can be given as 


TL=2(20logr+« 7) 


Equation (20) 


where o, =a, + O,, G, IS given by equation (17) and q, is the attenuation coefficient due to the 
presence of scatterers. The factor of 2 arises as the sound energy must travel from the transducer 
to the "target" and return to the receiver. Scattering losses due to the presence of bubbles and 


sediment will be covered in a later section. 


D. RESONANT BUBBLES 
Air bubbles suspended in water are subject to resonation due to excitation of the bubble by 
an incident sound wave. The frequency of bubble resonance is dependent upon the radius of the 


bubble according to the relationship 





ge ee 

he police 
2a a) 

Equation (21) 
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where f, = resonant frequency of the bubble in Hz, y = ratio of specific heats of the gas in the bubble 
(a value of 1.4 for air), p = density of water and P, = hydrostatic pressure. For air, at a depth of d 


feet under the surface, Urick [1983] simplifies this equation to 


2528 1 +0.03d 


Equation (22) 


As stated before, and by Urick [1983], only bubbles near resonant size will have an effect on the 
speed of sound within the water column. At the ultrasonic frequencies utilized by the CASP system, 
1.32 MHz and 5.2 MHz, the resonant bubble size, as determined using equation (22) above, is found 
to be approximately 2.3 um and 0.627 um, respectively, for near surface measurements. These 
bubble radii are far below the size range which has been resolved in previous field measurements, 
and even if bubbles do exist at the resonant size, viscous effects will probably prevent resonance. 
Viscous damping is most prevalent at small resonant radii (1 - 10 pm), but becomes less dominant 


at radii above 100 pum (Leighton [1994]). 


E. SCATTERING EFFECTS BY OCEAN BUBBLES 

Target Strength is defined by Urick [1983] as the echo intensity of incident acoustic energy 
measured one meter form an underwater target. In essence, target strength, TS, is a measure of the 
acoustically reflective characteristics of a target. In this case, the targets of interest are clouds of 
bubbles in the water column. For bubbles at resonance, the scattering cross section, o,, and the 
extinction cross section, 0, are approximately three orders of magnitude greater that the geometric 
cross section, (Leighton [1994]). This can be advantageous when using acoustic vice optical 


measurement methods. Leighton defines the extinction cross section, o , as the loss of energy from 


50 


an acoustic wave resulting from the presence of bubbles. It is given by the relationship 


4na’sd 
ka 


a 
(fF if? -1)? +8? 


Equation (23) 


For acoustically "large" bubbles, where a >> 1, small bubble theory does not apply. In this region, 


bubbles have an acoustic cross section equal to their geometric cross section, (Urick [1983]) 


Oo =1Ta 
Equation (24) 


Quite obviously, the volume scattering strength term is dependent upon the size of the scatterer, in 
this case, the radius of the bubble. If the bubble population is not made up of uniform sized bubbles, 
the scattering cross section, o, must be determined through integration over the range of bubble radii 
within the cloud 


6 ,=] n(a)o(a)da 


Equation (25) 


where n(a) is the density distribution of bubbles in units of bubbles per cubic meter per micron of 


radius, a , o(a) is the radius dependent scattering cross section and da is the integration interval. 


F. VOLUME REVERBERATION 
Volume reverberation is defined by Urick [1983] as the intensity of sound energy returned 


to the receiver when a sound source is projected into a volume containing scatterers. These 


SI 


scatterers can take the form of suspended sediment, zooplankton, phytoplankton or, in this case, air 
bubbles. The greater the number of scatterers, the greater the amplitude of backscattered energy. 
Volume reverberation is dependent not only upon the transmission loss through the medium, and 
the scattering strength of the scatterers within the volume as noted before, but also upon the 


ensonified volume 


Equation (26) 


where c 1s the speed of sound within the medium, t 1s the pulse duration of the transducer, D 1s the 
half power beam width of the transducer and r is the range to the ensonified volume. Table (5), 


below shows the applicable parameter of the Coherent Acoustic Sediment Flux Probe system. 


ae Transducer Wave- Farfield 
Radius, cm length, cm Range, cm 
degrees 
5202 |o3i7__| 0.028 


Table (5). CASP Monostatic Acoustic Beam Characteristics 
(after Anderson [1994]) 
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The above values are computed for a sound speed, c = 1480 m/s, in fresh water and a pulse duration, 
t = 2.26 x 10° seconds, (Anderson [1994]). Table (6), after Kinsler, et al [1982] shows values of 


sound speed for fresh and saltwater, for comparison. 


Se 







Temperature, degrees C Speed, c, m/s 
<a o. M48 
1500 


Table (6). Sound Speeds in Fresh and Salt Water 
(after Kinsler, et al [1982}) 


G. BACKSCATTERED ENERGY 


Urick [1983] defines volume reverberation level, RL, as reverberation power measured at the 


hydrophone. Reverberation is system specific and is given by the relationship 


RL = SL - TL,, - TL,, + SS + 10logV 
Equation (27) 


where SL is the source level of the transducer, TL,, and TL, are the transmission loss terms for 
source to scatterer and scatterer to receiver, respectively. For a monostatic sonar system, these two 
terms are identical. Although the CASP system can operate in either a monostatic or bistatic 
configuration, the sound ray paths have equal length, so the transmission loss terms can be 
combined into a single term. SS is the scattering strength of the scatterers present within the 
ensonified volume V. 

Echo level, EL, is defined as the intensity of the echo reflected from the target measured in 


the water at the hydrophone. It is given by the relationship 


EL =SL-2TL+ TS- DT 
Equation (28) 


Urick defines the combination of equations (27) and (28), EL - RL, as the echo to reverberation 
ratio. It is the difference between the strength of the echo and the strength of the reverberation due 
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to scatterers. It is independent of both the source and the medium in which the sound 1s transmitted. 
Combining all terms discussed above, the amount of backscattered energy received by the 


transducer head of the CASP system, EL - RL, becomes 
EL-KL=BS ,=1S—(SS+10l0gV) 
Equation (29) 


where BS, is defined as the sound energy backscattered to the receiver by a given volume containing 
scatterers such as bubbles. This value is still dependent upon the equipment in use, however, the 
ensonified volume term, shown in equation (26), is dependent upon the bin size, source transmission 
level, and the sensitivity and directivity of the receiver package. Target strength for a single bubble 
is given by the relationship 


Oo 
TS=10log— 
t 


Equation (30) 


If the scattering cross section o, 1s integrated over the entire distribution of bubbles, as noted 
previously for the volume scattering term, S,, a target strength can be obtained for an entire cloud 


of bubbles. 


H. VOLUME BACKSCATTER MODEL 

In order to obtain a theoretical comparison for laboratory and field measurements, a 
computer based acoustic model was written using MATLAB. Published bubble density 
relationships were used in order to determine a theoretical value for backscattered energy, and void 
fraction measurements. For this model, a variety of bubble densities, n(a), were used. These were 
presumed to be reasonable representations of the expected surf zone size distribution and spectral 


density. Several bubble size spectra have been reviewed in Section II.B.1. Medwin and Breitz 
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[1989] found that for bubbles below 60 jum in radius n(a) ~ a~, and for bubbles of radius greater than 
60 um, n(a)~ a”’. Averaged over the entire distribution, Medwin and Breitz found that n(a) = 7.8 
x 10°a”’. N(a) has the units of bubbles per cubic meter per micron radius. Cartmill and Su [1993] 
found that the number of bubbles per cubic meter per micron, n(a) ~ 2% while Su, Ling and Cartmill 
[1988] found the dependence to be n(a) ~ a“°. Figure (20) depicts bubble density from 10 to 500 
um using Medwin and Breitz's average density relationship. Figure (21) depicts a plot of theoretical 
backscattered energy values expected during laboratory experiments using the CASP system's 5.3 
MHz transducer, over a range of zero to 1.2 meters, neglecting the backscattered energy received 
as a result of the test vessel’s presence within the large tank. Figure (22) shows the backscattered 
energy predicted for a bubble cloud contained within the test vessel. For these calculations, acoustic 
system parameters are shown in Table (5) for the 5.3 MHz CASP frequency. Environmental and 


bubble size spectral parameters are shown below in Table (7). 


Sound Speed 1500 m/s 
Bubble Spectral Density n(a) = 7.8 x 10° (a *’) 


Table (7). Volume Backscatter Model Parameters 






The radius increment, da, shown in Table (7) above, is used in equation (25) in order to determine 
the scattering cross section of the modelled bubble population. Figure (23) shows a theoretical 
attenuation curve for the modelled bubble population. For this model, the theoretical void fraction 
was found to be 8.05 x 10°. This value corresponds well with published void fraction values. 
Leighton [1994], from Monohan, et al [1987] stated that estimates of void fraction under breaking 
waves suggest values of 10% % to 10° %. Lamarre and Melville [1992] from Walsh and Mulhearn 
[1987] and Farmer and Vagle [1989] found void fraction measurements to be within the range of 
10° to 10°. It should be noted that Monohan, Walsh and Mulheam, and Farmer and Vagle 
conducted open ocean measurements, rather than surf zone measurements. Table (8) shows a 


summary of field measurement conditions for comparison purposes. 
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Lamarre and Melville 4.5 - 14.9 m/s above 0.3 % 


Table (8). Field Void Fraction Measurement Conditions 






Lamarre and Melville's impedance probe was incapable of measuring void fractions lower than 
0.3%. The data shown above is a clear illustration that void fraction decreases rapidly with depth. 
At a depth of 10 meters, entrained bubbles due to wind action are likely to be very small, as larger 
bubbles would have risen to the surface before deep entrainment was accomplished. Near the 
breaking surface of the wave, as in Lamarre and Melville's experiment, the void fraction is expected 
to be much higher, with a broader spectrum of bubble sizes. Of interest, Lamarre and Melville 
[1992] state that at void fractions above 0.3%, temperature effects on the impedance of the 


Surrounding water are negligible. 
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Figure (20). Modelled Bubble Spectral Density, 10 to 500 um Radius 
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Figure (21). Modelled Backscattered Energy for 5.3 MHz CASP Frequency 
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Theoretical Backscatter Level, 5.3 MHz 
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Figure (22). Modelled Backscattered Energy for 5.3 MHz CASP Frequency in Test Vessel 
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Figure (23). Estimated Attenuation Curve 
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V. INITIAL EXPERIMENTAL TESTING 


A. SINGLE NITROGEN BUBBLE TESTS 

In an effort to visualize the laser light scattering effects of air bubbles in water and 
to realize potential illumination problems, a simple, controlled experiment was conducted. 
The apparatus used included the laser and mirror arrangement previously discussed in 
addition to the camera and microscope assembly used to videotape the results. Bubbles were 
produced using a single hypodermic needle fed by a pressurized nitrogen tank. The needle 
was submerged 1n the test volume with the exit point positioned just above the vertex of the 
two mirrors suspended below the laser and was placed at an angle in order to ensure that the 
needle itself would not interfere with either the laser light fan or the free rise of the bubbles. 
The flow of nitrogen from the tank was controlled using a metering valve to reduce the flow 
so that single bubbles were viewed 1n each frame. With a turbulent cloud of air bubbles in 
the test volume, it was very difficult to distinguish individual bubbles and therefore difficult 
to determine the spectral density of the population. With a single bubble, the specular 
reflection pattern could be easily studied. Additionally, the hypodermic needle produced 
a stream of uniformly sized bubbles which simplified the study of the reflective patterns. 


Equipment parameters are shown in Table (9). 


Metering Valve Inlet Bubble Flow Rate Camera Shutter Speed 


Pressure 


(1000 see 


Table (9). Single Nitrogen Bubble Test Equipment Parameters 





Figure (13) showed the orientation of the laser head in relation to the mirrors. The laser was 


fitted with a beam spreader which produced a fan of light 0.3 mm deep spread 15 degrees 
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either side of centerline. This arrangement essentially functioned as three light sources, with 
the laser illuminating bubbles from above and the two mirrors illuminating the bubble from 
below. Theoretically, this would produce three "hot spots” on the surface of each bubble. 
The beam spreader significantly reduced the laser light intensity, which required selection 
of a shutter speed which would provide sufficient illumination of the bubbles. It was found 
that shutter speeds slower than 1/1000 th of a second resulted in smearing of the images of 
bubbles as the bubbles rose. 

Resultant images of bubbles showed a distinct pattern of three bright dots on the 
surface of the bubble. These "triples" represented the reflected images of the laser light 
incident upon the surface of the bubble. The image was quite similar to that shown in Figure 
(10), from Figure (4) of Su, et al [1994]. In their image, three dots, produced by 
photographic strobes, appeared on surface of the bubble. The spacing of these dots was then 
used to determine the size of the bubble. 

With a single bubble present within a viewing frame, it is simple to perform the 
required measurements. Due to a bubble's position within the light fan, some would only 
display one or two "hot spots", rather than three. Obviously, a single hot spot could not be 
used to determine the bubble's size as detailed by Su, et al [1994]. Two "hot spots” could 
be used in this manner to provide a somewhat less robust result. However, with a cloud of 
bubbles present within a viewing frame, there was no simple method to separate individual 
bubble triples either visually or through the use of computer automated image analysis. In 
cases where an individual bubble displayed more than one "hot spots", the spots were often 
of different intensities. This was found to be a function of the onentation of the laser 
apparatus in relation to the camera. In order to achieve uniform illumination, the camera 
viewing plane and light fan had to be precisely parallel with the fan of light produced by the 


laser. To overcome these limitations a different approach had to be undertaken. 
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B. SINGLE REFLECTION METHODS 

There was no guarantee that each individual bubble would display three "hot spots" 
of equal intensity due to its position in the light beam. Statistically, this would pose a 
problem if a small number of bubbles were imaged and a large number of those displayed 
only one or two “hot spots". In single bubble tests, this theory was easily tested. Visual 
analysis of captured video images revealed that even in this controlled bubble flow, the 
number of hotspots displayed varied from bubble to bubble. For this reason, and given the 
turbulent nature of the test vessel during stirring, a single source illumination method was 
selected for use with a 90 degree scattering angle. 

During the laboratory runs described above, the brightness of the captured bubble 
image was dependent upon both the intensity of the incident laser light, and the shutter speed 
of the camera. In order to increase the resolution of a single "hot spot", to reduce image 
smearing, the shutter speed of the camera was reduced to 1/4000 th second. To further 
increase the resolution, the magnification of the microscope was doubled, which required 
a quadrupling of the incident light intensity to sufficiently illuminate the surface of the 
bubble. In order to increase the intensity of laser light, the beam spreader was removed to 
produce a narrow beam which was of much higher intensity than the fan of light produced 
with the beam spreader installed. Although the light energy distribution across the beam 
was considered to be Gaussian, this unmodified beam had a nominal diameter of 0.35 mm, 
and was considered to be cylindrical. The increased light intensity produced by this beam 
allowed for the use of a 1/4000 th second shutter speed selection on the camera. This would 
prove critical during bubble illumination in a turbulent tank in order to reduce video image 
smearing caused by illumination of high velocity bubbles. 

Figure (24) shows a plan view of the relative orientation of the equipment used 
during final experimental runs, which will be discussed in greater detail in Chapter VI. As 
can be clearly seen, there was very little unused space within the test vessel during 


experimentation. This was a primary factor in selecting a ninety degree scattering angle 
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rather than any other. The camera / microscope assembly was onented so that the viewing 


axis was perpendicular to the axis of the laser beam inside of the vessel. This onentation 


took into account the refraction of the beam as it passed through both the plexiglass walls 


of the tank and vessel and the water. 
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Figure (24). Laboratory Apparatus Tank Orientation 
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C. IMAGE PROCESSING 
The illuminated bubble cloud was imaged with a microscope and electronically 


shuttered video camera. Specifications for both the camera and microscope are shown in 


Tables (10) and (11). 


Cohu 4915 Monochrome HAD to 1/10000 sec 


Table (10). Video Camera Specifications 













Infivar Video Scope 1: 1.067 


Table (11). Microscope Specifications 


The aspect ratio is based on the viewing plane dimensions (512 x 480 pixels). Initially, a 
commercial image processing software package was sought to analyze video images in order 
to size and count the bubbles and to estimate a size spectrum. Several different packages 
were investigated, but each of these packages were rejected due to limitations in the 
scripting language needed to automatically size the bubbles. The software was required to 
automatically size and count a wide range of bubble radii within the test vessel, often sizing 
several objects within the same frame simultaneously. In order to achieve the required 
flexibility, a C-callable routine was developed "in-house" which wrote object sizes to an 
ASCII file which could later be analyzed using a MATLAB algorithm. 

Recorded video images were played back through a Matrox MVP framegrabbing 
package to capture individual video frames. Captured frames were analyzed for the 
presence of bright areas within the frame above a specified threshold level. These recorded 
images were measured and counted using a C-callable algorithm which measured both the 


width and height of the images within a captured frame. Careful selection of a threshold 
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level was required in order to ensure an accurate depiction of the spectrum of bubble 
images. Low threshold values resulted in a great deal of "noise" present within the frame 
while too high a value resulted in very few detected objects. Secondary external reflections, 
internal refractions and reflections are ignored using this method as the intensities of these 
bright areas are significantly lower than the primary hotspots. This is shown quantitatively 
in Table (4), from Van de Hulst [1957]. 

The sizing routine sampled roughly 100 frames per minute, freezing each frame and 
performing the searching and sizing algorithm. While the sizing routine measured both 
width and height, it was found that the height was an unreliable value due to the presence 
of blank horizontal lines of video data interlaced between each valid line of data when the 
camera was used in the shuttered mode. As described in Section I.B.8, the laser illuminated 
bright areas within each frame were not images of the entire exterior surface of the bubble, 
but were the specular reflection of the laser light off the exterior surface of the bubbles as 


seen by the video camera and microscope assembly. 
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VI. EXPERIMENTAL RESULTS 


A. BUBBLE PRODUCTION MEDIUM 

Experimental runs were conducted with both fresh and seawater in the test vessel. 
As different bubble spectra were expected, both fresh water and salt water runs were 
conducted. The plan was to evaluate data sets to not only confirm the produced bubble 
spectrum, but to determine the acoustic effects of bubble clouds within the test vessel. 
Using the chosen bubble production method of above water jet entrainment, bubbles were 
illuminated using a 45 milliwatt, 532 nanometer laser and counted and sized with the 
algorithm discussed in Chapter V. Optical Backscatterance Sensors provided a qualitative 
means to analyze the resident bubble population. Void fraction was estimated through the 
use of a conductivity cell. Additionally, the Coherent Acoustic Sediment Flux Probe 
(CASP) was used to produce acoustic backscatter levels for each bubble population. A 


detailed description of each of the above noted experimental steps follows. 


1. Fresh Water 

Fresh water experiments used fresh deionized and filtered water. Following 
deionization, the water was filtered through a 20 um and a | pm filter, in succession. These 
precautions were taken in an effort to minimize the amount of contaminants present within 
the test vessel. These contaminants could have taken the form of suspended particles, biotic 
particles, oils and other surfactants which would each have a potential effect on the 
produced bubble spectrum, and the backscatter properties of the medium. Additionally, the 


laser imaging system could not discriminate between bubbles and suspended particles. 


2. Simulated Sea Water 
Salt water runs were conducted utilizing simulated sea water. Salt water runs were 
necessary in order to simulate the conditions experienced in the surf zone as closely as 


possible. Sea water would effect the conductivity of the medium, as well as the distnbution 
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and size of the resident bubble population, as discussed in Chapter [V. Bubbles produced 
in salt water are, in general, more numerous and more persistent than in fresh water. Table 
(12), below, shows the constituent breakdown of the sea salt used to create the simulated sea 
water. The salt used in the experimental process was produced by the Lake Products 


Company of Maryland Heights, Missouri, under the product name, "Sea Salt". 


MgCl, - 6H,O 26.46 
















SrCl, - 6H,O 0.095 
a ee 0.007 


Table (12). Chemical Constituents of Sea Salt 


According to the manufacturer, "Sea Salt" 1s a simulated sea salt mix containing elements 
found in natural sea salt in quantities greater than 0.004%. In order to produce a simulated 
sea water, 41.953 grams of sea salt were added to each liter of fresh water to produce sea 
water, or for the test vessel, 600 grams of salt were added to approximately 14.5 liters of 
deionized fresh water. Table (13) shows the conductivity readings obtained using an 


AutoSal machine to determine the salinity of the water used during simulated sea water runs. 


68 





ee 
Table (13). AutoSal Salinity Test Results of Vessel Water 







23.08 


It should be noted that all three above trials were conducted at 24 degrees centigrade. 
Although the salinity of the sample tested was much lower than actual seawater, it was 
sufficient for the conductivity cell to operate effectively. As stated previously, deionized 
water and the sea salt were used in an effort to reduce the amount of contaminants within 
the test vessel, rather than actual sea water which could introduce biologic and other 


particles into the test volume. 


B. BUBBLE POPULATION VARIATION 

Two experimental runs were conducted using salt water, and one using fresh water. 
Fresh water was expected to produce larger bubbles when compared to the finer, smaller 
bubbles produced in salt water. In order to effectively validate the results, two salt water 
runs were necessary. During each run, four bubble production rates were attempted by using 
a needle valve to control the flow of water from the pump through the above water jet 
apparatus. Water flow was controlled in one half turn steps. Following approximately ten 
minutes zero bubble unstirred data, the vessel was stirred for five minutes prior to bubble 
production. The valve was opened one half turn, producing a low flow rate through the jet. 
With this jet flow, the corresponding bubble population produced within the test vessel was 
allowed to stabilize for ten minutes. The valve was then opened an additional half turn for 
an additional ten minutes. This procedure was continued for two more half turn steps, each 
for ten minutes. It was hoped that each of these steps would produce a bubble cloud of 


increasing population density without altering the spectrum produced. A discussion of the 
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spectrum produced follows in section E. Table (14) shows point indicators and the 
corresponding needle valve setting. These indicators apply to all time series plots in 


following sections for saltwater experimental runs conducted on 03 and 07 August 95. 





Table (14). Time Series Point Indicators 


C. OPTICAL BACKSCATTER RESULTS 

Optical backscatter results were obtained using the Optical Backscatterance System 
(OBS) manufactured by the D & A Instrument Company of Port Townsend, WA. A 
specification sheet for this system is included as appendix (A). The OBS system is designed 
to be used for turbidity measurement in water and is not advertised as having the capability 
to detect bubbles. Bubbles, apparently, do not have the proper opacity to be effectively 
detected by the system. Despite this, it was felt that this system would provide another 
means of detecting changes in bubble population density, even if this was only on a 
qualitative basis. It was not expected that this system would provide accurate quantitative 
data with respect to bubble populations produced within the test vessel. Figures (25) and 
(26) show the optical backscatter voltage time series produced during two separate salt water 
experiments. The upper plot is the time series decimated to a 1 hz sampling rate, while the 


lower one shows a time series decimated to a 0.05 Hz sampling rate. This has the effect of 
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removing perturbations caused by large bubble transients. In each of these fi gures, the open 


circles denote times during the run where the jet flow was increased. 


OBS ch 30 data, 1 Hz sampling 


2 O 
ro 
> 
Ss 
a 
6 -0.5 
¢p) 
Mm 
O 





> 
0 oe 0.4 0.6 0.8 1 12 1.4 


OBS ch 30 data, 0.05 Hz sampling 
0.1 


© 


OBS Output, Volts 





0 0.2 0.4 0.6 0.8 1 1.2 1.4 
Elasped Time of Run, hrs 


Figure (25). Optical Backscatter, 03 August 95 


TI 


OBS ch 30 data, 1 Hz sampling 


OBS Output, Volts 





0 0.2 0.4 0.6 0.8 1 lez 1.4 


OBS ch 30 data, 0.05 Hz sampling 


0.2 
L E 
S 0.1 A B 6 D 
—o O O O O O 
ra 
£ 0 
O 
8 
0 -0.1 

0.2 

Or 0.4 0.6 0.8 1 1.2 


Elasped Time of Run, hrs 


Figure (26). Optical Backscatter, 07 August 95 
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As can be clearly seen from these two figures, there is a definite increase in the signal 
following the initial opening of the valve, approximately 0.2 hours into the experimental 
run in Figure (25) and approximately 0.25 hours into the run in Figure (26). An increase is 
again apparent at point (B). Subsequent increases in flow rates, at points (C) and (D) are 
less distinct in these figures. This not only indicates that the OBS system is at least partially 
capable of detecting bubbles within a volume, but more importantly, that the bubble density 
was, in fact, changing with increases 1n flow rate through the jet, at least through point (C). 
The large perturbation in Figure (26), near point (E) was caused by high speed agitation of 
the test volume with the laboratory stirrer mechanism. This was done in an effort to check 
the effects of a completely saturated environment on the sensors in use. During this 
agitation, the vessel contained a very high void fraction and was nearly opaque. The degree 
of the opacity is apparent from the OBS signal. 


D. CONDUCTIVITY CELL RESULTS 

Void fraction was estimated from the conductivity of the fluid within the test vessel 
during each experiment. Since the presence of air in the fluid decreases the electrical 
conductivity as the number of bubbles within the vessel increases, the fluid's conductivity 
decreases. An inductive conductivity cell produced by the Falmouth Scientific Company 


was used 1n the tank. 


1. Cahbration 

The conductivity cell was calibrated with an autocell to understand the relationship 
between counts (representing digitized voltages within the CASP), volts, and mmohs / cm, 
units of conductivity. A specification sheet on the conductivity cell is included as appendix 
(B). Calibration data for the conductivity cell used in the experiments is shown below in 
Table (15), and voltage calibration of the analog input was determined by connecting a 


precision voltage source and measuring the resulting digitized counts. Each voltage setting 
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Conductivity Calibration Resmits 67a 35 


Instrument type: 
Model: FSI 
Sensor :Conductivity 

Sensor interface: 1637800592 


Conduetivi ty 


Standard serial number: Y2 

Calibrator: BRUCE 

Comments: PRE COPE EXPERIMENT 

Point # X 
Volts 


1.10590 
1.1058@ 
-48980 
. 48980 
03170 
.03160 
-.49822 
-.49830 
= 29170 
Se dhice 


Ow oan nnarawnn- 


— 


Conductivity Calibration Results 8/8/95 


Instrument type: 
Model: FSI 
Sensor: Conductivity 

Sensor interface: 1637800592 


Conductivity 


Standard serial number: Y2 
Calibrator: BRUCE 
Comments: PRE COPE EXPERIMENT 


Linear model: Y=A+B+#x 
A= 35.0041021701 
B= 15.0240280901 


X in Volts Y in deg C (IPTS68) 


Si 


Y 
Conductivity, mmohs/cm 


.68652 
.69000 
- 46390 
.48950 
-22100 
sy ONS) 
ome 
.37460 
.45112 
-42490 


X 


Se eo oqce ddd 9 & 


Manufacturer: 
Serial number: 


Sensor serial 


Bridge serial 


stddev 


.Q0E+@0 
. 8Q@E+02 
.Q@0E+02 
. Q0E+00 
.Q@GE +00 
.Q0E+200 
.Q0E+92 
.90E+2@2 
.Q0E+O0 
.Q@GE+00@ 


Manufacturer: 
Serial number: 


Sensor serial 


Bridge serial 


NPS 


1494 


number : 


number: iQ 


Seo eooeseooo sd 9 


NPS 
1494 


number: 


number : 


stddev 


. Q0E+202 
. GVE+2Q 
. 9G@E+O0 
.QGE+00 
.Q0E+200 
. 2@GE+2@ 
. QVE+Q@ 
. ZVE+Q0O 
. Q@0E+0Q 
.Q9GE+O0 


1494 


1027 


Quality 


Seo ge og gogo oso D9 


Fit parameter R=(explained sum of squares/total sum of squares?) ##1/2 = 


s893869754126 


Table (15). Conductivity Cell Calibration Results 
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was maintained for 15 seconds before selecting the next step. With a known input voltage, 
a corresponding response was determined in an effort to obtain a conversion factor between 


voltage and counts. The results of this calibration are shown below in Table (16). 





Input Volts, DC Cell Response, counts 
43,9883 


Po -12830 6393 


] 
Z 
3 
4 
5 









6385.18 


Table (16). Conductivity Cell Voltage Calibration Data 


Taking the initial offset of -44 counts at an input of 0 volts into account, the relationship 
between input voltage and cell response tn counts is linear. Calibration testing resulted in 


the following relationship 


C=( 46.74mmohs/cm \ volts Neonne) 


volts 6383 .2counts 
Equation (31) 


where C is the conductivity of the fluid. As a second test of the conductivity cell's response 
within the test vessel, several readings were taken with the conductivity cell in a variety of 
positions, both in and out of the water. Results of this test are shown in Table (17). It 
should be noted that the readings were fairly consistent throughout the test. Of note, the 


proximity of the edge of the cell to the wall or other obstructions had little or no effect on 
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the conductivity reading. However, if the open end of the cell was placed near another piece 
of equipment or the vessel wall, there was a marked effect, quite similar in magnitude to the 


reading obtained during the high turbulence stirring portion of the test. 





a | 
Te 












Position / Condition 


Table (17). Conductivity Cell Position Test Results 






Notes - * Side of cell touching wall, ** Opening of cell touching wall 


2. Conductivity Results 
The conductivity cell was used to estimate total void fraction data for each of the 


produced bubble populations. Void fraction is determined from the relationship 





Equation (32) 


where o,, is the conductivity of the water, o., is the effective conductivity of the water 
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containing bubbles and a is the void fraction of air within the sample volume. 
determination of the conductivity and effective conductivity of the water will be discussed 
in the next section. Figures (27) and (28) show conductivity cell signals plotted against 
elapsed experimental run time. Once again, the open circles denote times when the flow 
rate of the entrainment jet was increased. Unlike the OBS data shown in Figures (25) and 
(26), there is no distinct decrease in the conductivity signal following flow increases. Figure 
(27) shows the most distinct change in the conductive nature of the medium with an increase 
in the bubble population after the second increase in flow rate, at point (B). Following that, 
there 1s a gradual decrease in conductivity, rather than distinct steps. It is likely that the 
conductivity cell in use experienced an unknown malfunction as the cell can resolve 0.002 
mmohs / cm changes in conductivity. This is confirmed by the fact that after bubble 
production had ceased, approximately one hour into each run, the conductivity cell readings 
did not return to the level registered at the beginning of the run. Initially, it was thought that 
this was a result of bubbles becoming trapped within the conductivity cell's center. This was 
discounted when the laboratory stirring apparatus was run at high speed momentarily. The 
result of this effort can be clearly seen in Figure (28), at approximately one hour. The 
conductivity reading immediately returned to the value registered at the end of bubble 
production. Data shown in Figure (28) is further corrupted by the presence of suspended 
particles which will be discussed briefly in the next section. Unfortunately, this prevented 


void fraction from being accurately estimated using conductivity cell data. 
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Figure (27). Conductivity Cell Results, 03 August 95 
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Figure (28). Conductivity Cell Results, 07 August 95 
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E. BUBBLE PRODUCTION AND SPECTRA RESULTS 

To confirm that the bubble clouds in the test tank represented oceanic spectra, a large 
effort was made to determine the bubble spectra in the test vessel. The method for sizing 
bubbles was based on an accurate measurement of the size of the specular "hotspot" on the 
bubbles surface (see Section If.B.8). From Baldy and Bourguel [1985], the relationship of 
this specular "hotspot" and the size of the bubble is 


eae = 0.06 A subble 
Equation (33) 


A major limitation of this method was the finite resolution imposed by the camera pixel size. 
A calibration of the camera and microscope was made by imaging a target within the test 


vessel, with results shown in Table (18). 


Target Width Frames Imaged Mean Target Width 
333 pixel 16.5 um / pixel 


Table (18). Video Calibration Results 












Using the relationship as given in equation (29), this would result in resolution of air bubbles 
only down to a radius of approximately 138 tum. It was initially hoped that this method 
would resolve bubble sizes down to 10 tum in radius. The lower limit of bubble radius was 
determined by the size of a single pixel, and the upper limit was restricted only by the 
nominal diameter of the laser beam. This resulted in an upper limit of over 2900 um in 
radius, with bubbles of larger radius being measured as having a 2900 pm radius. 
Regardless of these restrictions, however, a portion of the bubble size spectrum should still 
have been resolved using this method, allowing comparisons of bubble spectra produced 


within the test vessel with measured ocean spectra reviewed 1n Sections II.B.1 and 2. 
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Figure (29) shows results of fresh water bubble production, while Figures (30) 
through (33) show results of salt water bubble production. Table (19) shows experimental 
results used in these plots. Bubble size spectral density, n(a), (bubbles / radius / meter’), is 


shown as a function of bubble radius. Spectral density is determined by the relationship 


n(a)= 





datVxi 


Equation (34) 
where N is the number of bubbles of a given radius, da is the radius increment, V is the 
sample volume and / is the number of frames of data which were processed. These plots 
also depict results published by Cartmill and Su [1993] for comparison. The most notable 
difference between these data sets is that our laboratory spectral densities are larger by over 
two orders of magnitude. This suggests that the bubble population in our experiments are 
far higher than the open ocean measurements reported in the literature. Additionally, the 
very small sample volume (4.247 x 10° m’) limits the minimum number of objects of a 
given radius illuminated during an experimental run. Ifa single bubble is detected within 
the sample volume during an experimental run, the bubble population automatically contains 
10,000 bubbles m®, when the sample is normalized to a cubic meter. This is clearly shown 
in Figure (30) which shows a single data point at 1400 microns radius. This plot shows the 
particle size spectrum present within the vessel prior to bubble production. It is likely that 
a single suspended particle was present, was illuminated and counted by the sizing routine 
and became represented as 10° particles upon normalization. The measurements do, 
however, show that the observed bubble size spectral levels are much higher than the open 
ocean values measured by Cartmill and Su [1993], and yet visually better represent bubble 


densities lighter than those observed in the surf zone. 


81 














Table (19). Video Bubble Image Data 
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# Bubbles/radius/meter’3 


Bubble Spectrum, Nozzle Setting 2, 16 Jun 95, Fresh Water 


lab data 
cartmill & su - fresh 
cartmill & su - salt 





Bubble radii, microns 


Figure (29). Freshwater Bubble Spectrum 


# Bubbles/radius/meter43 





Bubble Spectrum, Nozzle Setting 0, 03 Aug 95, Salt Water 


lab data 
cartmill & su - fresh 
cartmill & su - salt 


Bubble radii, microns 


Figure (30). Saltwater Bubble Spectrum, 03 August, No Bubble Production 


84 


# Bubbles/radius/meter’3 


O 
Je 


Bubbie Spectrum, Nozzle Setting 0.5, 03 Aug 95, Salt Water 


lab data 
cartmill & su - fresh 
cartmill & su - salt 





Bubble radii, microns 


Figure (31). Saltwater Bubble Spectrum, 03 August, Minimum Flow 
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# Bubbles/radius/meter’3 


Bubble Spectrum, Nozzle Setting 2, 03 Aug 95, Salt Water 


lab data 
cartmill & su - fresh 
cartmill & su - salt 





Bubble radii, microns 


Figure (32). Saltwater Bubble Spectrum, Maximum Flow 
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# Bubbles/radius/meter43 


Bubble Spectrum, Nozzle Setting 0, 07 Aug 95, Salt Water 


lab data 
cartmill & su - fresh 
cartmill & su - salt 





Bubble radii, microns 


Figure (33). Saltwater Bubble Spectrum, 07 August, Still Conditions 
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The most disturbing aspect of the laboratory measurements shown in Figures (29), 
(31) and (32) is that the slope of the curve in all runs is much steeper than that of Cartmill 
and Su. This is probably an artifact of the measurement method. Despite the choice of an 
HAD, CCD camera, extremely bright pixels could "bleed" over into adjacent pixels, making 
the imaged hotspot larger than it is. This phenomenon will be more prevalent in smaller 
bubbles than in larger ones, and would have the effect of steepening the curve in the region 
of smaller radii, producing the shape shown. Corruption of spectral data is also apparent in 
Figure (33) which shows still water conditions with no bubble production and quite a few 
suspended particles within the vessel. These particles, which could have been biotics 
produced in the tank during a weekend, could have had an adverse effect on conductivity 


readings and corresponding void fraction estimation. 


F. ACOUSTIC DATA 


Coherent Acoustic Sediment Flux Probe data was used to determine the effects of 
the bubble cloud within the test vessel upon a transmitted acoustic signal. As discussed in 
Chapter IV, the scattering strength of the bubbles within the test vessel 1s a function of their 
geometric cross section integrated over the population's distribution. For this reason, it was 
expected that the strength of the backscattered energy would be greater with more bubbles 
in the vessel, and that the attenuation of the CASP signal through the cloud would increase. 
This was borne out through analysis of the data obtained during the two salt water 


experiments. 


1. Backscattered Energy 


The amplitude of the acoustic energy received by the CASP, which is backscattered 
by the bubble population resident within the test vessel, 1s dependent upon the scattering 
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strength (effective cross section) of the bubble population. It was initially hoped that the 
optical measurement methods employed during these experiments coupled with the bubble 
sizing algorithm would provide a quantitative basis for determining the effects of the 
produced bubble size spectra on the backscattered energy. Due to the fact that the optical 
measurement method and conductivity cell did not function as well as planned, no 
quantitative verification of the size spectrum was available. As a result, backscattered 
energy could only be analyzed in terms of the backscatter level, attenuation at two 
frequencies and the optical backscatter levels without explicit knowledge of the test bubble 
size spectrum in an effort to determine an attenuation coefficient within the vessel. Figures 
(34). and (35) show backscattered energy levels at the front and back of the vessel, 
respectively, for 03 August experimental runs, for nozzle settings A through E. Figures (36) 
and (37) show similar data for 07 August experimental runs. Energy levels recorded at the 
front of the vessel are of higher magnitude and display much more distinct changes when 
bubble production rates are increased. At back wall of the vessel, CASP signal strength 1s 
attenuated by the bubble cloud so that raw backscattered levels are nearly constant, 
regardless of the bubble population. Attenuation results will be discussed in the next 


section. 


2. Attenuation 


Acoustic attenuation effects were discussed in Chapter IV. Attenuation coefficients, 


a, for each separate bubble population produced was determined as a function of range 
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within the test vessel. Acoustic pressure as a function of range, P(r), 1s given by the equation 


-(a +a )Ar 


: 
P(r)=P (—)e 
r 


Equation (35) 
where P,, is the reference acoustic pressure taken at the first bin inside the test vessel. The 
reference range, r,, 1s also taken as the first bin inside of the vessel, while the range, r, runs 
from the front to the back of the vessel. For both of these values, the bins containing the 
walls of the vessel were not used in the calculation. The attenuation coefficient for water, 
Q,, 1S a predetermined quantity based on the density of the water and frequency in use. The 


attenuation coefficient due to the presence of scatterers, a, 1s given by the relationship 


aw =-(In(—+)+a,) 
Pr 


Equation (36) 


A least squares fit of these data points produced a line whose slope is the attenuation 


coefficient for that particular population. Table (20) shows a summary of the attenuation 


values obtained. 





Table (20). Attenuation Coefficients 
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This data clearly shows that for runs 3 and 4, there is very little difference in the bubble 
population resident within the test vessel. Companson with the attenuation curve obtained 
using the bulk backscattering model in Chapter IV shows a reasonably good agreement given 
the bubble spectral density used in the model (shown in Table (7)). 


3. Suspended Sediment Attenuation Comparison 


Kohanowich [1995], during studies of acoustic attenuation due to suspended 
sediment, confirmed that attenuation of the CASP signal at 1.3 MHz was linearly dependent 
upon the concentration of scatterers present within the test vessel according to the following 


relationships 


= 
a =-0.0857—-— 
(g/l) 


Equation (37) 


and the received acoustic power at the 1.3 MHz transceiver was related to the concentration 


of sediment in grams / liter (up to a concentration of 725 gm/1) 1s 


5.2 x 10% (g/l) / counts? 
Equation (38) 


Therefore, it is a simple matter to determine a sediment concentration which would be 


required to produce an attenuation coefficient equivalent to those measured due to the 


presence of bubbles. Results of these calculations are shown in Table (21). 


=) 


Date / Point Indicator Attenuation Coefficient, Counts” Equivalent Sediment 
a, (x 10°) Concentration, (g/l) 


0.913 





Table (21). Suspended Sediment Equivalence 


As can be seen from the above data, even high bubble concentrations, as were produced 
during experimental runs 3 and 4, are equivalent to relatively low concentrations of 
sediment. In open ocean measurements, Lamarre and Melville [1992] found void fractions 
near the surface of breaking waves to be greater than 0.3 %. In the surf zone, void fraction 
near the surface of breakers could reasonably be expected to be much higher. This value 
would attenuate rapidly with depth, while the concentration of entrained sediment would 
increase to a maximum near the bottom. In a position within the water column where both 
bubbles and sediment are present, the majority of backscatter received by the CASP would 
be from entrained sediment. Figure (38) shows two plots, both of which compare the 
attenuation effects of bubbles and sediment. In each plot, *'s denote laboratory 
measurements, the solid line is a linear least squares fit of this data, and the dotted line show 


attenuation information obtained using equations (33) and (34) for entrained sediment. 


2 


Clearly, the difference in the slopes indicate that a discriminator could be developed to 
determine the difference in backscatter levels due to the presence of both air bubble and 


sediment scatterers with the surf zone. 
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Figure (34). Backscattered Energy, Front of Vessel, Bin 14, 03 August 
Seawater, Nozzle Settings A- E 
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Figure (35). Backscattered Energy, Back of Vessel, Bin 20, 03 August 
Seawater, Nozzle Settings A - E 
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Figure (36). Backscattered Energy, Front of Vessel, Bin 15, 07 August 
Seawater, Nozzle Settings A - E 


96 





1.4 


log10 Counts*2 


Wo 


o> 
On 


Oo) 


Sd 
on 


on 


i 
on 


oS) 


Bin20 Signal for CASP beam 1 


A B C D E 


G2 0.4 0.6 0.8 1 1.2 
Elapsed Run Time, Hrs 


Figure (37). Backscattered Energy, Back of Vessel, Bin 20, 07 August 
Seawater, Nozzle Settings A- E 
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Vil. CONCLUSIONS 


The initial goal was to develop a means by which a bubble population within a test 
vessel could be accurately counted and sized in order to produce a spectral density similar 
to that which could be reasonable expected to be encountered in the surf zone, although 
there were no surf zone observations with which to compare. This representative population 
would then be utilized to determine the ultrasonic backscatter characteristics of different 
bubble densities. This data would allow further determinations to be made concerning a 
method to discriminate between the signal produced by acoustic energy scattered from 
suspended sediment and the signal scattered by entrained bubbles. The development of this 
method of measurement was hampered by several influences. Time constraints, complexity 
and availability of measurement apparatus, and the physical restriction imposed by 
laboratory equipment all played significant roles in not only the selection of equipment to 
be used, but also the eventual end of the experimental process. 

A significant effort was made to utilize both optical measurement methods and void 
fraction measurements to assess the bubble population produced within the test vessel. At 
very low bubble densities, optical measurement methods are very accurate and the data is 
simple to analyze. As demonstrated during single bubble production testing, measurement 
and counting of the produced bubbles using laser illumination and video image analysis was 
quite easy. Bubbles of the size range which are expected to be encountered in the surf zone 
effect the conductivity of the water they are suspended in only when taken as a bulk 
measurement. Individually, they are not large enough to have a significant effect on the 
conductivity. At high bubble densities, regardless of the bubble spectra, conductivity cells 
function well, and changes in spectral density can be detected. High spectral densities, 
however, have the effect of making the water nearly opaque, rendering optical measurement 
methods ineffective. In order to use both of these methods together in the same volume, a 


"middle ground" spectral density had to be produced that would allow a low enough bubble 
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density for the use of laser optics to measure the individual bubbles, yet high enough for 
void fraction to be determined using conductivity. This was nearly accomplished. 

The chosen method of bubble production, above water jet entrainment, produced, at 
least visually, a reasonable bubble spectrum. In view of the shortcomings of the method of 
measurement, however, this spectrum could not be verified quantitatively. 

The method of bubble measurement used appeared to have the ability to measure and 
size individual bubbles accurately, however, it proved to have one major drawback when 
used to determine spectral density. Geometrical optical approximations are a well accepted 
means of determining the size of air bubbles in water, but the extremely small sample 
volume presented by the laser and camera apparatus used limited the smallest number of 
objects measured during and experimental run. Additionally, the potential “smearing” of 
pixels as presented in video images of illuminated bubbles may have induced errors in the 
measured size of the "hotspots" produced on very small bubbles. Although this method was 
only capable of sizing bubbles down to 140 um, this was not viewed as a limitation. 
Theoretically, this method should still have been capable of resolving some portion of 
published spectral measurements which were used as a benchmark for this expenment. 
Limitations imposed by the small sample volume and the resolution capabilities of the video 
equipment prevented this. 

Qualitative analysis of void fraction results showed that although the bubble density 
was apparently altered through control of the flow through the jet, the equipment in use 
suffered some sort of malfunction which prevented accurately detecting subtle changes in 
void fraction. Calibration testing of the conductivity cell revealed very little difference in 
cell response for zero bubble conditions and saturation of the vessel with bubbles. This is 
regarded as significant in that variations in spectral density throughout expenmental runs 
were visually quite different, but resulted in imperceptible changes in the conductivity of the 
medium as registered by the conductivity cell. Additional difficulties may have ansen due 


to the lack of a completely homogeneous bubble cloud and the restrictive physical 
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conditions of the vessel which prevented orienting the cell in a manner which would have 
allowed for free flow of bubbles through the center. 

Acoustic backscatter results revealed nothing unusual. Studies of scattering strength 
and attenuation due to clouds of bubbles are well documented. This series of experiments 
produced results which did nothing more than validate theories and studies published by 
several investigators. One interesting result, however, was the validation that the Coherent 
Acoustic Sediment Flux Probe was quite capable of detecting relatively subtle changes in 
bubble density. the CASP proved to be far more sensitive than any of the other sensors used 
in these experiments. 

Despite a significant effort, it is clear that the methods chosen for bubble 
measurement and counting were not the best methods available. From the outset, it was felt 
that an Acoustic Resonator Array would be the simplest, most accurate and least invasive 
of previously published methods available. Unfortunately, efforts to obtain or manufacture 
a resonator proved to be too difficult. In addition to this, the size of the test vessel precluded 
the use of an array large enough to have a sufficient measurement range of bubble radii. 
Further efforts in this area should include the use of an acoustic resonator in concert with 


a conductivity cell. 
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APPENDIX A. SPECIFICATION SHEET FOR OPTICAL 
BACKSCATTERANCE SENSORS 
OBS-3 & OBS-3B 


SUSPENDED SOLIDS & MECHANICAL DATA 
TURBIDITY MONITOR OBS-3 


FEATURES 
71 (181) 


OBS-3 monitors have field-proven Optical BackScatterance Cans 
(OBS®) sensors, low-power analog circuitry, and switch- 
selectable gains with fine adjustment. These self-contained 
sensors are ideal components for CTD’s, water quality 
monitoring and process contro! systems and can monitor 
suspended matter in the clearest to most turbid waters in 
rivers, lakes, and the ocean. 


12 (31) 








The Model OBS-3 includes an underwater housing rated to 1 (25) 45 (114) nT y7() ee 
2,000 or 5,000 meters and ts available with an optional 4 mA 
to 20 mA two-wire transmitter. A cable assembly is required to 
connect the Model OBS-3 to a power supply/data logger. Con- 
sult the manufacturer about your cable requirements. 


inches tmull 





The Model OBS-3B is designed for the system intergrator. It in- OBS-3B 
cludes a threaded sensor with flexible, color-coded #30 AWG SENSOR 
leads for connections to the circuit module. Color-coded #30 

AWG leads are also provided for connecting the circuit 

module to an OEM power supply/data logger. 


APPLICATIONS 


Vv OEM component in oceanographic, water quality, 
and process-control instrumentation 

Water and wastewater monitoring systems 
Laboratory measurements 

River gauging 

Vv Sediment transport studies 


SPECIFICATIONS 


DYNAMIC RANGES 






1.2 (31) 


qd<4< 


0562 +002/ - 000 









o-ring surfaces 


4 
} 
l 
1 
4 
J 
4 


Mnoidity .......2.- _. 002 to 2,000 FTU' 

UC) (a ae .... 0.1 to 5,000 ppm (me/I) 

S20) ears .... 2 to 100,000 ppm (me/l) ge ON 
FREQUENCY RESPONSE ........ 10 Hz 10 (254) 

AMBIENT LIGHT REJECTION _.. . Optical filtering and to logger/power supply 


synchronous detection 


TEMPERATURE COMPENSATION ... Solid state temperature 
transducer 


39 (98) 


| 338 (8539) 









NONLINEARITY2 





6). 
Turbidity (formazin, 0 to 2000 FTU) . 2.0% © 
Mud (0 to 4,000 me/l) ..... me 205 8 
Sand (0 to 60,000 mg/l) .... ..... 3.5% 2 
DRIFT 
ine 22... < 3.5% in Ist 2,000 hrs CIRCUIT tf 051 113) 
Temperature .... 2... .. 2.2... 0.05% /°C = 
SUPPLY VOLTAGE COMPLIANCE .. 250 pV/V MODULE 12020 en 
SETTLING TIMES ‘©5O6O O ue 
POWCIUD «2.229540. ales 
25 °C Step Change in Water 
Temperature ...... ee 15S 
OUTPUT SPAN ..... Pe Ot0 3 V 1 Coaxial Shield 4 Detector (Coaxia! lead) G Gain 
OUTPUT IMPEDANCE ee 200 Ohms 2 Common (green) 5 Temperature Sensor (orange) Z Offset 
RMS NOISE AT O FTU < 50 wW 3 IRED Anode (red) 6 JRED Cathode (brown) T Temperature Coe 
B 

POWER REQUIREMENTS Serre 7 1Onls ya ni inches Im. 
OUTPUT FILTER ... . . . 20 Hz (-3 dB) 
PHYSICAL DIMENSIONS DETER EEO ARO | 
Sensor . ~ length eS eran) TERMS: COD (except with prior credit approval) 

diameter 0.7 in (18 mm) 
Housing length 4.5 in (115 mm) INSTRUMENT 40-A Seton Road 

diameter . _. 1.2 in (31 mm) COMPANY Port Townsend, WA 98368 
Weight . a oe 04 Ib (0.17 kg) - ea Phone: 800-437-8352 

submerged . 0.16 tb (0.07 kg) 103 5 tee Cee ee 385-0272 ‘ 

7206 -0460 

WORKING DEPTH... .. 2,000 mor 5000 m aa a 


1 Formazin Turbidity Uni 
2 Maximum deviation of response from a straight line. expressed as a 
percentage of the calibration range. OBS sensors are protected by US Patent 4.841.157 





APPENDIX B. SPECIFICATION SHEET FOR CONDUCTIVITY CELLS 


SPECIFICATIONS MEASUREMENT: 





CONDUCTIVITY: | RANGE 0 -> 64 mmho/em ( 7 
ACCURACY +/- .025 mmho/com (.0025 S/m) 
STABILITY +/- .005 mmho/cm/month (.0005 S/m) 


SE 
SS 




















TEMPERATURE: | RANGE -5 -> 35 Celsius (ITS-90) | 
RECURReY Met) _050 Celsius 


| | STABILITY +/- .002 Celsius/Month | 


SPECIFICATIONS ELECTRICAL & MECHANICAL: 





= —_ 


PROBES : Pa lnicue ite cHeonirhac Inductivereonductay ity | 
Sensor, Platinum Temperature Sensor 





SIZE SENSOR: Max. O.C. 45 mm LENGTH 163 mm 
SIZE CARD: Width 50.8 mm, Length 101.6 mm, Height 20 mm 


MOUNTING: 3/2°SR SS PPR Witworth Thread) 


O-RING: 2-123 
MATERIAL: Sensor (active area): Ceramic 


OPERATING: 1000 PSIA (650 dBAR) ABSOLUTE MAXIMUM 
PRESSURE: 


S vpc +/- 10% JUMPER J1 INSTALLED 
P=) 4 eve JUMPER J1 REMOVED 
LOAD CURRENT @ 5 mAMPS 


100 mSeconds from power up 
SAMPLING: 100 mSeconds Between Channel Selections | 
Tis CONS TA! je 10 mSecond Conductivity, 100 Seconds Temp. 


CONNECTOR: Molex 6 Position Housing Type #2695 









INPUT OCGrC : tar” (oD Care aoe ype max’ S "S VDC | 


OUTPUT : OM >) 2 eseVvolEeS=Nominal ** | 


Mameawcile EaAccory EGr eperonal output voltage ranges. 
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SENSOR OUTLINE DRAWINGS 


2-123 O-RING 


4 8 
3/4 BRITISH STANDARD PIPE THD 


CERAMIC /DELRIN COND SENSOR 14 THOS/INCH 





STANDARD C-T SENSOR 


1.80 OO x .91 10 xX 1.80 LC 


EXPANSION BLADOER 
(DEEP RATED) 





LARGE CONDUCTIVITY ONLY (EXTERNAL TEMP) SENSOR 
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